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EXECUTIVE  SUMMARY 


OBJECTIVES 

Demonstrate  the  value  of  high-performance  visualization  techniques  for  the  method- 
of-moment  modeling  of  electromagnetic  field  radiation  and  scattering  for  antennas  in  a 
complex  environment.  Improve  the  utility  of  computational  techniques  used  in  electro¬ 
magnetic  analysis,  with  special  emphasis  on  ship  antenna  design. 

APPROACH 

Method-of-moments  is  one  of  the  most  commonly  used  techniques  in  computational 
electromagnetics.  It  is  especially  useful  for  many  of  the  frequencies  of  interest  in  the 
electromagnetic  topside  design  of  Navy  ships.  The  method-of-moments  technique  involves 
solving  a  set  of  coupled  integral  equations  to  determine  the  currents  on  a  structure. 

Computational  electromagnetics  can  be  thought  of  as  a  three-step  process:  problem 
definition,  computation,  and  solution  description.  The  drawback  to  using  a  computational 
technique  such  as  method-of-moments  is  in  the  effort  required  at  each  of  the  three  steps. 
Input  models  must  be  extensively  validated  before  calculations  are  performed.  Enormous 
quantities  of  output  data  are  generated,  including  currents  on  all  wire  segments,  near-field 
contours,  and  far  fields.  Advanced  visualization  techniques  can  assist  in  input  validation 
and  the  rapid  interpretation  of  output  data. 

The  computational  electromagnetics  code  used  was  the  Numerical  Electromagnetics 
Code  -  Method  of  Moments  (NEC-MoM).  Electromagnetic  visualization  was  performed  on 
a  Silicon  Graphics  Incorporated  (SGI)  4D/320GTXB  workstation. 

CONCLUSIONS 

Visualization  tools  now  exist  for  performing  analysis  and  validation  at  all  three  steps 
of  a  method-of-moments  design  procedure.  These  tools  have  greatly  improved  the  utility 
of  computational  techniques  used  in  electromagnetic  analysis  and  are  now  being  used  to 
support  topside  ship  design  projects  for  the  Navy. 
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1.0  INTRODUCTION 


This  report  presents  the  findings  of  a  Navy-funded  Independent  Exploratory  Develop¬ 
ment  (lED)  investigation.  This  investigation  applied  advanced  visualization  techniques  to 
an  existing  computational  electromagnetic  (CEM)  code.  The  goal  was  to  demonstrate  that 
high-performance  visualization  can  improve  the  utility  of  computational  techniques  used 
in  ship  electromagnetic  designs. 

The  CEM  code  used  was  the  Numerical  Electromagnetics  Code  -  Method  of  Moments 
(NEC-MoM).  The  NEC-MoM  code  was  not  run  on  the  visualization  workstation.  Instead,  it 
was  run  on  a  high-performance  computing  (HPC)  platform,  a  Convex  C-220,  with  only 
the  NEC-MoM  input  and  selected  output  files  ported  to  a  Silicon  Graphics,  Incorporated 
(SGI)  machine  for  pre-  and  postprocessing  visualization. 

The  visualization  platform  used  was  an  SGI  4D/320GTXB,  representing  the  most 
sophisticated  visualization  hardware  available.  Two  visualization  codes  were  developed  as 
a  result  of  this  project;  viewjiec  and  viewjnom.  The  programming  was  done  using  the  C 
language  calling  SGI  Graphics  Library™  (GL)  functions. 

This  report  is  organized  along  the  following  lines:  Chapter  2  describes  computational 
electromagnetics,  the  method-of-moments  technique,  and  the  need  for  visualization. 
Chapter  3  discusses  techniques  for  visualizing  data.  Chapter  4  discusses  the  approach 
used  for  this  lED  project,  including  hardware  and  software  options.  Chapters  5,  6,  7,  and 
8  give  detailed  descriptions  of  the  visualization  products  that  were  developed  as  a  result 
of  this  project.  Chapters  9  and  10  provide  observations  and  conclusions. 

2.0  COMPUTATIONAL  ELECTROMAGNETICS 

Computational  electromagnetics  is  that  branch  of  electromagnetics  that  routinely 
involves  using  a  computer  to  obtain  results.  It  is  a  complementary  tool  to  the  classical 
techniques  of  experimental  observation  and  mathematical  analysis. 

The  CEM  application  of  interest  to  the  authors  of  this  report  is  naval  ship  electromag¬ 
netic  design.  The  specific  area  of  interest  is  antenna  design.  Key  technical  parameters  for 
antenna  performance  include  impedance,  near  fields,  coupling,  and  far-field  radiation 
patterns. 

2.1  ELECTROMAGNETIC  MODEL 

The  electromagnetic  model  uses  transfer  functions  derived  from  Maxwell’s  equations 
(figure  1).  The  inputs  to  the  transfer  functions  include  an  environmental  description  of  the 
problem  as  well  as  the  specified  input.  The  environmental  description  is  defined  by  both 
the  electrical  and  the  geometrical  properties  of  the  structures  and  the  space  in  which  they 
reside.  The  input  is  the  specified  excitation,  which  can  be  either  a  voltage  source  applied 
to  an  antenna  or  a  plane  wave  impeding  on  the  defined  structure.  The  outputs  from  the 
transfer  functions  are  the  induced  currents  on  the  structures  or  the  impedance  of  the 
antenna.  These  induced  currents  can  be  used  to  determine  both  the  near  and  far  fields. 
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For  naval  ship  electromagnetic  ..'sign  the  near  fields  are  of  interest  in  determining 
coupling  between  antennas  as  well  as  hazards  to  personnel,  fuel,  and  ordnance.  Far  fields 
are  used  to  determine  the  performance  of  the  electromagnetic  system. 


ENVIRONMENT  DESCRIPTION 
ELECTRICAL,  GEOMETRICAL 


INPUT 


(SPECIFIED 

EXCITATION) 


TRANSFER 

FUNCTION 


DERIVED  FROM 


MAXWELL’S 

EQUATION 
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(INDUCED  SOURCES: 
-NEAR  FIELDS 
-FAR  FIELDS) 


Figure  1.  Electromagnetic  model. 


2.2  NEC-MoM 

The  numerical  electromagnetic  code  was  developed  for  antenna  modeling.  NEC 
includes  NEC-MoM,  NEC-basic  scattering,  and  NEC-reflector  antenna  codes  (Li,  Logan, 
and  Rockway,  1988).  All  these  codes  have  been  validated  and  extensively  documented. 
NEC-MoM  is  a  development  effort  of  the  Lawrence  Livermore  National  Laboratory 
(Burke  and  Pogio,  1981;  Burke,  1989). 

To  model  complex  structures,  NEC-MoM  uses  an  electric-field  integral  equation 
(EFIE)  for  wires  and  a  magnetic-field  Integral  equation  (MFlB)  for  surfaces.  The  EFIE  is 
well  suited  for  thin  wires,  while  the  MFIE  is  attractive  for  modeling  smooth  surfaces  that 
bound  large  closed  volumes.  The  EFIE  can  also  be  used  to  model  surfaces,  by  using  a 
wire  grid  to  represent  a  surface.  The  EFDE  wire  grid  is  preferred  for  surfaces  that  are  not 
smooth  and  that  enclose  volumes  that  are  not  large.  For  the  task  described  in  this  report, 
only  the  EFIE  technique  was  used. 

The  thin-wire  approximation  is  applied  to  the  EFIE  to  reduce  the  equation  to  a  scalar 
integral  equation.  Transverse  currents  and  transverse  variation  of  the  axial  current  on  the 
wire  are  neglected.  The  boundary  condition  on  an  electric  field  is  enforced  in  the  axial 
direction  only.  The  approximations  made  in  NEC-MoM  are  valid  as  long  as  the  wire 
radius  is  much  smaller  than  the  wavelength  and  the  wire  length.  The  result  is  a  set  of 
coupled  integral  equations  for  the  currents.  As  a  practical  matter,  these  approximations 
limit  the  application  of  NEC-MoM  to  resonance  and  below  for  a  given  structure. 
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The  integral  equations  are  solved  numerically  by  the  collocation  form  of  the  method  of 
moments,  which  involves  expanding  the  unknown  current  in  a  summation  of  basis  func¬ 
tions  and  enforcing  the  equality  of  weighted  integrals  of  the  fields.  This  reduces  the  inte¬ 
gral  equations  to  a  matrix  equation  (Harrington,  1968).  The  weighting  functions  are 
chosen  to  be  delta  functions.  This  collocation  method  results  in  a  point  sampling  of  the 
fields.  Wires  are  divided  into  short  segments,  with  a  sample  point  at  the  center  of  each 
segment. 

The  matrix  equation  that  results  from  the  moment  method  is  solved  by  LU  decomposi¬ 
tion  (Press  et  al.,  1986).  This  determines  the  currents  on  the  wire  segments.  From  the 
currents,  all  electromagnetic  parameters  of  interest  can  be  derived:  impedance,  coupling, 
and  near  and  far  fields. 

2.3  NEED  FOR  VISUALIZATION 

Computational  electromagnetics  can  be  thought  of  as  a  three-step  process:  problem 
definition,  computation,  and  solution  description.  The  drawback  to  using  a  computational 
electromagnetics  code  such  as  NEC-MoM  is  in  the  effort  required  at  each  of  the  three 
steps.  Preparing  the  input  model  for  NEC  and  evaluating  the  output  can  be  an  overwhelm¬ 
ing  task.  Performing  calculations  can  be  exceedingly  time  consuming  for  all  but  the  most 
simple  structures. 

To  conserve  computing  resources,  input  models  must  be  extensively  validated  before 
calculations  are  performed.  Enormous  quantities  of  output  data  are  generated,  including 
currents  on  all  the  wire  segments,  near-field  contours,  and  far  fields.  One  of  the  goals  of 
this  lED  project  was  to  demonstrate  how  advanced  visualization  techniques  could  be  used 
to  assist  in  input  validation  and  the  rapid  interpretation  of  output  data. 

3.0  VISUALIZATION  TECHNIQUES 

During  this  lED  project,  decisions  had  to  be  made  regarding  how  to  effectively  visual¬ 
ize  the  data  of  interest.  In  many  cases,  there  were  several  options.  It  was  often  difficult  to 
decide  which  display  method  would  have  the  most  utility  for  the  final  user.  The  products 
of  this  project  are  now  being  used  in  support  of  ship  EM  design  projects.  This  application 
of  these  products  will  quantify  utility. 

3.1  COLOR  CODING 

One  of  the  main  problems  encountered  was  in  deciding  how  to  assign  color  coding. 
Color  coding  is  an  art  unto  itself.  It  was  often  difficult  to  decide  whether  color  coding 
should  be  done  using  a  continuous  range  of  colors  or  using  a  discrete  set  of  colors. 

There  are  three  components  to  color.  The  RGB  color  model  interprets  these  compo¬ 
nents  as  the  three  colors  (red,  green,  blue)  used  by  the  CRT.  By  varying  the  amount  of 
each  color  component,  the  full  range  of  screen  colors  can  be  achieved.  The  SGI  machine 
has  24-bit  color,  with  8  bits  for  each  component.  Each  component  is  an  integer  value 


3 


be  .vten  0  and  255.  In  this  manner,  the  console  is  capable  of  displaying  over  16  million 
different  colors.  Of  course,  it  is  doubtful  that  the  human  eye  can  distinguish  all  these. 

Humans  do  not  perceive  colors  in  the  same  manner  that  the  console  displays  them.  To 
humans,  it  is  more  intuitive  to  use  the  HSV  color  model  (hue,  saturation,  and  value).  The 
HSV  model  is  based  on  the  intuitive  appeal  of  the  artist’s  tint,  shade,  and  tone.  Hue  has  a 
value  between  0  and  360.  Saturation  and  value  range  between  0  and  1.  Algorithms  exist 
for  translating  between  the  different  color  models  (Foley  et  al.,  1990). 

At  one  point  during  this  ED  project,  it  was  proposed  that  various  data  components 
could  be  encoded  into  the  different  components  of  the  HSV  color  model.  For  example, 
for  complex  data,  the  magnitude  could  be  encoded  in  the  hue  and  the  phase  could  be 
encoded  in  the  value.  This  idea  was  found  to  be  not  feasible  for  two  reasons.  First,  the 
user  became  -verwhelmed  by  the  amount  of  information  contained  in  slight  variations  in 
displayed  color.  Second,  the  color  printer  was  totally  incapable  of  reproducing  anything 
but  major  variations  in  color,  so  one  could  not  obtain  a  hard-copy  output  of  the  visualiza¬ 
tion  display. 

For  most  data  sets,  a  discrete  color-key  coding  system  was  found  to  provide  the  most 
useful  visualization  of  the  data.  Seven  bins  were  chosen  and  the  data  were  linearly 
assigned  to  the  bins.  The  color-key  assignments  were  changed  several  times.  A  final  deci¬ 
sion  was  made  to  use  colors  that  gave  the  best  contrast  when  printed  out  on  the  color 
printer.  Phase  is  displayed  as  a  continuous  range  of  hue,  with  saturation  and  value  set  to 
1.0. 

3.2  COMPLEX  DATA 

In  most  cases,  we  found  that  trying  to  display  both  components  (real  and  imaginary  or 
magnitude  and  phase)  of  a  complex  data  set  in  one  image  confused  the  user  unnecessar¬ 
ily.  A  decision  was  made  to  allow  the  user  to  display  complex  data  in  side-by-side  win¬ 
dows  if  both  components  needed  to  be  viewed  simultaneously. 

3.3  VECTOR  DATA 

Current  data  and  field  data  are  vectors.  A  straightforward  method  of  displaying  a 
vector’s  orientation  was  never  developed.  The  current  vector  was  defined  relative  to  the 
wire’s  direction.  It  was  felt  that  the  direction  of  the  current  did  not  have  as  much  value  as 
the  magnitude  and  phase  of  the  current. 

4.0  APPROACH 

For  advanced  visualization,  decisions  had  to  be  made  regarding  hardware  and  soft¬ 
ware.  The  authors  were  fortunate  to  have  access  to  an  SGI  4D/320GTXB.  This  machine 
represents  the  most  sophisticated  visualization  hardware  currently  available.  Several  visu¬ 
alization  packages  were  available  on  the  SGI  machine,  including  PV-Wave,  apE,  Explorer, 
and  the  SGI  Graphics  Library.  The  available  high-  level  languages  included  C  and 
Fortran. 
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4.1  SGI  4D/320 


The  SGI  4D/320GTXB  is  one  of  the  SGI  POWER  Series^^  line  of  computers.  It  can 
be  configured  with  up  to  eight  CPUs;  the  one  used  for  this  project  has  two  33-MHz  CPUs. 
It  has  the  following  CPU  performance  ratings  (SGI,  1991):  59  MIPs  (VAX  Dhrystone 
MIPS),  20  MFLOPS  (DP  Linpack  1000x1000),  41  SPECmarks.  The  GTXB  has  48  bits 
color  and  24  bits  Z  buffer.  Its  graphics  performance  is  rated  as  400K  vectors/sec,  150K 
triangles/sec,  and  lOOK  polygons/sec.  The  system  used  for  this  project  is  configured  with 
64-MB  memory  and  a  19-inch  console  monitor. 

4.2  SOFTWARE  PACKAGES 

Several  graphics  software  packages  were  available  on  the  SGI  computer.  Due  to  time 
and  monetary  constraints,  no  investigation  was  made  of  purchasing  other  packages.  A 
decision  was  made  to  work  with  the  available  packages.  The  available  packages  were: 
PV-Wave,  apE,  Explorer,  and  the  SGI  GL. 

PV-Wave  is  a  very  powerful  visualization  package.  Unfortunately,  it  did  not  seem  to 
run  properly  on  the  SGI;  moreover,  it  was  only  suited  to  graphically  displaying  data  sets. 
It  had  no  3D  model  rendering  capabilities  that  would  be  needed  for  the  proposed  BED 
project. 

apE  is  a  software  toolkit  for  visualization.  It  was  developed  by  the  Ohio  Superco¬ 
mputer  Graphics  Project  starting  in  1987.  Scientific  and  engineering  data  can  be  proc¬ 
essed  and  viewed  as  plots,  color  images,  and  three-dimensional  objects.  Unfortunately,  it 
was  cumbersome  to  use.  3D  objects  had  to  be  described  using  the  apE  format.  It  was  very 
difficult  to  translate  from  a  NEC  input  file  to  apE.  For  this  reason,  apE  was  not  used  for 
this  lED  project. 

Explorer  was  developed  by  SGI  to  fully  use  the  extraordinary  capabilities  of  the  SGI 
visualization  computers.  It  is  similar  in  design  to  apE,  but  it  is  much  less  cumbersome  to 
use.  Explorer  is  a  system  for  creating  powerful  visualization  maps,  each  of  which  com¬ 
prises  a  series  of  small  software  tools,  called  modules.  A  map  is  a  collection  of  modules 
that  carries  out  a  series  of  related  operations  on  a  data  set  and  produces  a  visual  repre¬ 
sentation  of  the  result.  Explorer  has  a  data  conversion  utility  for  moving  data  between 
Explorer  and  other  data  formats.  A  module  builder  allows  custom  modules  to  be  created. 
Unfortunately,  Explorer  became  available  only  near  the  end  of  the  BED  project  year,  so  it 
was  only  used  for  a  small  portion  of  the  project.  There  was  also  concern  about  future 
portability,  since  Explorer  is  only  available  for  the  SGI  system. 

GL  is  a  set  of  graphics  and  utility  routines  that  provides  high-  and  low-level  support 
for  graphics.  The  routines  can  be  called  from  either  C  code  or  Fortran  code.  Though  quite 
primitive,  GL  allows  one  to  access  all  the  powerful  visualization  capabilities  of  the  SGIs, 
including  3D  drawing,  Gouraud  shading,  device  polling,  double  buffering,  coordinate 
transformations,  hidden  surface  removal  (z-buffering),  lighting,  pick  correlation,  and  tex¬ 
turing,  depending  on  the  hardware’s  capabilities.  SGI  provides  an  enormous  number  of 
demo  programs  that  can  be  easily  modified.  GL  provides  a  straightforward  way  to 
develop  the  capabilities  required  of  this  lED  project. 
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4.3  NEC  DATA  SETS 


The  NTEC  input  data  set  is  created  by  a  program  called  NEEDS  (Li,  Logan,  and  Rock¬ 
way,  1988).  The  data  set  is  an  ASCII  file  written  in  an  archaic  format  with  each  line 
representing  a  “card.”  The  first  two  characters  on  each  line  refer  to  an  alphabetic  code 
that  determines  what  the  data  on  that  line  represent.  For  example,  “CM”  is  a  comment 
line  and  “GW”  is  a  wire  description  line.  The  output  from  NEC  is  an  enormous  ASCII 
text  file  in  a  format  thrt  is  very  difficult  to  read  from  within  another  computer  program. 
For  that  reason,  a  pars  ir  program  is  needed  to  extract  the  data  of  interest  and  put  such 
data  in  a  software-readable  format.  A  preprocessing  filtering  program,  filter,  was  written 
in  Fortran  to  put  the  input  and  output  data  in  a  form  that  would  allow  easy  display  of  the 
data  of  interest.  Appendix  A  contains  a  source  code  listing  of  this  program. 

Filter  prompts  the  user  for  the  name  of  a  NEC  input  file.  The  user  must  enter  the  full 
filename.  Then  the  user  is  asked  whether  there  is  an  associated  currents  file.  If  the  answer 
is  yes,  the  user  is  prompted  to  enter  the  full  name  of  the  currents  file.  Finally,  the  user  is 
prompted  for  the  name  for  the  output  file.  This  file  must  have  the  extension  “.geo”  (for 
geometry) . 

Filter  expects  the  format  of  the  currents  file  to  be  as  follows;  two  lines  followed  by  the 
data  listed  as  “magnitude  phase”.  The  current  data  for  each  segment  are  on  their  own 
line.  The  last  line  in  the  file  must  be  “-1.234  -1.234”. 

5.0  VISUALIZATION  PROGRAM  STRUCTURE 

The  visualization  programs  that  were  developed  for  this  lED  project  were  written  in 
the  ANSI  C  language  using  the  SGI  Graphics  Library  to  access  the  visualization  tools.  The 
program  to  visualize  the  problem  definition  and  solution  description  is  named  viewjiec, 
and  a  source  code  listing  of  it  is  in  Appendix  B.  The  program  evolved  as  a  result  of 
feedback  from  users.  However,  the  basic  structure  using  multiple  windows,  pull-down 
menus,  and  transformations  was  determined  from  the  beginning. 

5.1  WINDOWS 

The  beauty  of  working  in  an  X-windows  type  environment,  such  as  that  on  the  SGI 
machines,  is  the  flexibility  it  affords  the  user.  That  flexibility  was  retained  in  the  develop¬ 
ment  of  viewjiec.  Each  data  component  is  displayed  using  3D  imagery  in  its  own  window. 
At  any  given  moment  during  the  running  of  the  program,  the  user  has  complete  control 
over  how  many  data  windows  are  displayed  as  well  as  their  sizes  and  locations.  In  addi¬ 
tion,  multiple  copies  of  viewjiec  can  be  launched  to  do  side-by-side  comparisons  of 
different  NEC  runs. 

5.2  MENUS 

Pressing  the  right-hand  mouse  button  while  the  mouse  cursor  is  within  any  viewjiec 
window  brings  up  a  “pull-down”  menu.  This  pull-down  menu  allows  the  user  to  select 
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which  data  windows  to  open  or  close,  which  transformation  is  active,  and  whether  a 
coordinate  axis  is  displayed.  Once  a  selection  has  been  made,  the  affected  windows  are 
updated. 

5.3  TRANSFORMATIONS 

Entire  books  have  been  written  on  the  display  of  3D  models.  This  section  will  just 
summarize  the  display  capabilities  available  in  view_nec.  The  available  transformations 
include  rotation,  translation,  and  zooming.  Each  3D  transformation  is  represented 
internally  in  the  computer  by  a  4  x  4  matrix.  The  IRIS  Geometry  Engines  transform  all 
geometric  data  (vertices  of  points,  lines,  and  polygons)  by  multiplying  each  vertex  by  the 
accumulated  matrices.  In  viewjiec,  the  transformations  are  actuated  with  the  left  mouse 
button  held  down  and  the  cursor  dragged  across  the  screen. 

There  are  three  basic  classes  of  transformations  that  can  be  carried  out  by  the  Graph¬ 
ics  Library;  projection  transformations,  viewing  transformations,  and  modeling  transfor¬ 
mations.  A  good  analogy  is  to  a  camera  with  a  versatile  lens.  Projection  transformations 
describe  the  type  of  lens  on  the  camera.  Viewing  transformations  determine  where  the 
camera  is  positioned  and  in  which  direction  it  is  pointed.  Finally,  modeling  transforma¬ 
tions  affect  the  location,  orientation,  and  size  of  the  3D  geometric  models  in  the  scene. 

There  is  often  more  than  one  way  to  carry  out  the  transformation  of  a  scene.  For 
instance,  instead  of  moving  the  camera  toward  the  object,  the  object  could  be  moved 
closer  to  the  camera.  However,  there  may  be  subtle  differences.  The  transformation 
methods  used  in  viewjiec  were  chosen  by  selecting  demo  programs  that  gave  the  desired 
effect. 

5.3.1  Zoom 

The  zoom  feature  allows  the  user  to  expand  a  selected  portion  of  the  model  for  a  more 
detailed  display. 

Zooming  is  achieved  internally  by  using  a  projection  transformation,  perspective. 
By  modifying  the  fovy  parameter,  the  field  of  view  in  the  y  dimension  is  modified.  This 
has  the  effect  of  making  the  scene  appear  closer  or  farther  away.  The  aspect  ratio 
between  the  field  of  view  in  x  and  the  field  of  view  in  y  and  the  distances  to  the  near  and 
far  clipping  planes  are  kept  constant. 

5.3.2  Rotate 

The  rotate  feature  allows  the  user  to  rapidly  change  the  orientation  of  the  model.  This 
permits  a  better  feeling  for  the  three-dimensional  nature  of  the  model  and  often  allows 
various  features  of  the  data  to  become  more  apparent. 

Rotation  of  the  scene  is  carried  out  internally  using  the  polarview  command. 
Polar  view  is  a  viewing  transformation.  This  command  moves  the  viewpoint.  The  dis¬ 
tance  from  the  origin  is  kept  constant  while  the  azimuthal  angle  in  the  x-y  plane  and  the 
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incidence  angle  in  the  y-z  plane  are  modified.  This  is  akin  to  moving  the  camera  over  a 
spherical  surface  surrounding  the  scene. 

5.3.3  Translate 

The  translate  feature  allows  the  user  to  move  the  model  up  or  dovm  or  left  or  right  in 
the  plane  of  the  screen.  This  shifts  the  origin  and  is  useful  before  zooming  in  on  a 
selected  portion  of  the  model. 

Translation  of  the  scene  is  carried  out  internally  by  using  the  modeling  transformation 
translate.  The  program  is  set  up  to  allow  only  translation  of  the  object  in  the  plane  of 
the  screen.  The  size  of  the  object  is  preserved. 

5.4  PICK  CORRELATION 

It  quickly  became  apparent  that  a  method  was  needed  to  selectively  choose  a  single 
element,  such  as  a  wire  segment,  from  the  3D  wire  object  display.  This  was  needed  for 
troubleshooting  as  well  as  linking  the  display  back  to  the  NEC  input  data  set.  The  middle 
mouse  button  was  chosen  for  this  purpose.  The  technique  used  is  known  as  pick  correla¬ 
tion.  The  Graphics  Library  provides  this  capability. 

Pick  correlation  identifies  objects  on  the  screen  that  appear  near  the  mouse  cursor. 
Information  about  these  objects  is  stored  in  a  buffer.  Viewjiec  uses  this  information  to 
allow  the  user  to  select  a  single  wire  segment  from  the  object  displayed  and  list  all  infor¬ 
mation  about  it.  A  sample  wire  parameters  window  is  shown  in  figure  2. 
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Figure  2.  Wire  parameters  window. 
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5.5  USER’S  GUIDE 


The  user  does  not  need  to  understand  the  internal  structure  of  viewjiec  that  is 
described  above  in  order  to  use  the  program.  Running  viewjiec  is  a  very  straightforward 
procedure. 

If  the  user  is  interested  in  viewing  only  the  geometry  description  products  or  wire 
currents,  then  just  the  file  with  extension  ''.geo”  is  required.  This  file  contains  the 
geometry  of  the  model  and  the  currents  on  the  wires  if  these  have  been  calculated  by  an 
NEC  run.  The  ''.geo”  file  is  generated  by  running  the  program  filter  as  described  above.  If 
the  user  is  interested  in  examining  the  near  or  far  fields,  then  additional  files  are  needed. 
The  far-field  data  are  contained  in  four  files:  *.etm  (magnitude  of  E-theta  component), 
*.etp  (phase  of  E-theta  component),  *.epm  (magnitude  of  E-phi  component),  and  *.epp 
(phase  of  E-phi  component).  The  near-field  data  are  contained  in  the  file  *  fid.  The  user 
must  generate  these  files  before  running  viewjiec.  For  the  required  format  for  each  of 
these  files,  consult  the  source  code  listing  of  viewjec  in  Appendix  B. 

To  run  viewjiec  type  "viewjiec”.  The  program  then  prompts  the  user  for  the  file 
name.  This  is  the  file  name  for  the  *.geo  file.  DO  NOT  TYPE  THE  .geo  EXTENSION! 
During  the  running  of  viewjiec,  the  program  will  automatically  access  any  additional  files 
needed  for  the  display  of  near  or  far  fields  if  requested.  All  the  files  must  have  the  same 
prefix.  ITie  extensions  must  be  assigned  as  outlined  above. 

During  the  running  of  viewjiec,  the  menu  can  be  brought  up  at  any  time  by  pressing 
the  right-hand  mouse  button  while  the  mouse  cursor  is  within  any  viewjiec  data  window. 
The  menu  allows  the  user  to  open  or  close  any  data  windows,  change  the  transformation 
mode,  or  toggle  the  coordinate  axis  display. 

The  transformation  mode  is  indicated  in  the  lower  left  comer  of  all  data  windows.  The 
transformation  is  carried  out  by  clicking  the  left  mouse  button  and  dragging  the  cursor 
within  any  data  window. 

Clicking  on  the  wire  model  with  the  middle  mouse  button  brings  up  the  wire  parame¬ 
ter  window  described  above.  The  currently  selected  wire  segment  flashes  until  it  is 
unselected  (by  either  clicking  the  middle  mouse  button  in  the  background  or  selecting  a 
new  wire  segment). 

Clicking  the  middle  mouse  button  in  the  near-fields  window  changes  the  near-field 
display  threshold  (this  is  described  below  in  the  solution  description  products  section). 

To  exit  viewjiec,  click  the  left  mouse  button  in  the  shaded  band  at  the  very  top  of  any 
display  window  and  select  "close”. 
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6.0  GEOMETRY  DESCRIPTION  PRODUCTS 

The  NEC  input  file  describes  the  modeled  object  (such  as  a  ship)  as  a  collection  of 
wires.  The  description  of  each  wire  includes:  a  tag  number  (for  identification  purposes), 
the  number  of  segments  on  the  wire,  the  (x,y,z)  coordinates  for  both  the  beginning  and 
end  points,  and  the  radius  of  the  wire.  Since  some  wires  are  tapered,  extra  information  is 
given  for  these  wires  to  describe  the  tapering. 

There  are  four  windows  that  can  be  opened  up  to  display  the  model’s  input  geometry. 
The  windows  each  display  one  of  the  following:  Wire  Segmentation  (segment  length  in 
meters).  Wire  Radius  (in  meters).  Segment  to  Radius  Ratio,  or  Wire  Connectivity  (none, 
one,  or  both  ends  connected).  All  of  these  windows  display  the  model  as  a  3D  wire  object. 
The  data  are  encoded  in  the  color  of  each  wire  segment  using  a  linear  color  assignment 
scheme.  A  color  key  is  displayed  in  the  lower  left  corner  of  the  window.  Figure  3  is  a 
sample  window  showing  wire  segmentation.  Without  color,  one  cannot  get  a  true  under¬ 
standing  of  the  “value”  of  the  display  to  the  user. 


a  SegmentuUon  •  Q 


1  h  t-  irUJUSl  KOTHTES:  I  <^30,  38q,  D.DD,  D.OO 


Figure  3.  viewjiec  wire  segmentation  display  window. 

7.0  SOLUTION  DESCRIPTION  PRODUCTS 

Solution  description  is  displayed  in  up  to  eight  windows.  These  eight  windows  include 
currents  on  the  wires  (real  component,  imaginary  component,  magnitude,  or  phase),  total 
near  field  (electric  field,  in  volts/meter),  z-component  of  near  field,  theta  component  of 
far  field,  and  phi  component  of  far  field. 
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7.1  CURRENTS 


Currents  on  the  wires  are  displayed  using  the  same  technique  as  was  used  for  the 
geometry  description  products:  current  is  color-coded  on  a  3D  wire  display  of  the  model. 
Since  currents  are  complex,  there  are  four  different  current  data  windows  that  can  be 
displayed:  real  component,  imaginary  component,  magnitude,  and  phase.  Most  users  have 
found  the  magnitude  window  to  be  most  useful. 

7.2  NEAR  FIELDS 

NEC-MoM  allows  the  user  to  calculate  the  near  field  at  selected  locations  surrounding 
the  model.  These  locations  are  usually  defined  by  a  3D  grid  surrounding  the  model.  The 
near  field  is  a  complex  vector.  Of  most  interest  is  the  z-component,  since  this  will  be  the 
component  that  affects  personnel  standing  on  the  surface  of  the  ship.  Several  techniques 
were  investigated  for  displaying  the  near-field  data  points. 

Viewjiec  displays  near  fields  using  a  “fog”  technique  in  which  the  density  of  activated 
pixels  in  the  image  is  linearly  proportional  to  the  field  intensity  at  the  nearest  calculation 
point.  Since  the  calculation  points  usually  form  a  3D  grid,  the  image  contains  square 
blocks  of  fog  varying  in  density  with  the  field  intensity.  The  fog  is  color-coded  in  a 
manner  similar  to  that  used  for  the  geometry  description  products.  The  3D  wire  model  is 
drawn  using  a  dark  gray  color  so  as  not  to  detract  from  the  near-field  display. 

The  near-field  windows  have  a  thresholding  capability.  By  clicking  in  the  window  with 
the  middle  mouse  button,  the  user  is  able  to  selectively  change  the  threshold  at  which  the 
fields  are  displayed.  This  allows  the  user  to  quickly  determine  what  areas  surrounding  the 
model  have  fields  above  a  particular  cutoff  level.  The  threshold  can  be  selected  to  be  any 
of  the  color-key  bin  levels,  HERP  (Hazardous  Electromagnetic  Radiation  for  Personnel), 
or  HERO  (Hazardous  Electromagnetic  Radiation  for  Ordnance).  HERP  and  HERO  are 
functions  of  the  frequency  at  which  the  NEC  model  was  run.  The  near-field  windows 
display  the  value  of  the  frequency.  HERP  and  HERO  values  are  shown  in  table  1. 

Table  1.  HERP  and  HERO  near-field  values. 


Frequency  Range,  MHz 

Maximum  E-Field,  V/m 

HERP 

0.01-3 

632.5 

3-30 

1897/(frequency) 

30-300 

63.25 

300-1500 

3.65  >$QRT  (frequency) 

1500-300,000 

141.4 

HERO 

0.1-1 

1-3.7 

100/f  frequency) 
100/(frequency^) 

3.7-10 

2 

10-1000 

0.5  X  (frequency®-®) 
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7.3  FAR  FIELDS 

There  are  two  windows  available  for  the  display  of  far  fields.  These  windows  allow  the 
user  to  display  either  the  theta  component  or  the  phi  component  of  the  far  electric  field. 

The  far  field  is  complex  data.  This  makes  display  of  it  somewhat  complicated.  A 
method  was  sought  of  displaying  both  the  magnitude  and  phase  simultaneously.  The 
method  selected  involves  displaying  the  fields  as  a  three-dimensional  surface.  The  dis¬ 
tance  from  a  point  on  the  surface  to  the  origin  is  proportional  to  the  field  magnitude  at 
that  point.  The  color  at  the  point  is  determined  by  the  phase  of  the  field  at  that  point. 
Gouraud  shading  is  used  between  points  to  allow  for  color  transition. 

8.0  MoM  MATRIX  DISPLAY 

There  has  been  some  speculation  that  visualizing  the  method  of  moments  matrix  that 
has  been  extracted  from  NEC  might  give  some  insight  into  the  wire  model’s  validity.  The 
MoM  matrix  is  a  complex  matrix  that  may  have  a  dimension  of  about  1000.  The  dimen¬ 
sion  of  the  matrix  is  equal  to  the  number  of  wire  segments  in  the  input  model.  This 
matrix  is  also  referred  to  as  the  impedance  matrix.  Two  methods  were  used  to  render  the 
data:  an  Explorer  map  and  a  Graphics  Library  program. 

8.1  EXPLORER  MAP 

An  Explorer  map  was  developed  to  display  the  magnitude  of  the  matrix.  This  map 
reads  in  the  data,  subsamples  and  crops  the  data,  extracts  the  magnitude  values,  and  then 
renders  the  data  as  a  three-dimensional  surface.  The  vertical  offset  of  a  point  on  the 
surface  (as  well  as  the  color  of  the  point)  is  proportional  to  the  magnitude  at  that  point. 

There  are  several  advantages  to  using  the  Explorer  map.  Explorer  allows  a  great 
amount  of  flexibility  in  assigning  color  mappings,  manipulating  the  rendered  3D  object, 
and  modifying  of  the  vertical  offset. 

There  are  also  several  serious  disadvantages  to  using  the  Explorer  map.  Foremost  is 
the  array  size  limitation  (64K),  which  does  not  allow  a  display  of  the  entire  matrix  for  a 
typical  ship.  A  typical  ship  matrix  can  have  as  many  as  a  million  elements.  The  array  can 
be  cropped  and  subsampled,  but  a  great  deal  of  data  are  lost.  This  makes  the  rendered 
graphic  almost  useless.  In  addition,  only  the  magnitude  of  the  data  can  be  displayed.  The 
phase  component  of  the  data  cannot  be  displayed.  Finally,  there  is  no  way  for  the  user  to 
determine  where  in  the  matrix  (i.e.,  row  and  column)  different  features  in  the  display  are 
located. 

8.2  GRAPHICS  LIBRARY  PROGRAM 

A  Graphics  Library  progreim,  viewjnom,  was  developed  to  allow  the  display  of  almost 
any  size  method-of-moments  matrix.  Its  structure  is  similar  to  that  of  viewjiec,  but  it  is 
quite  a  bit  simpler.  Viewjnom  allows  data  to  be  displayed  in  three  display  modes:  2D 
(surface),  ID  (wire),  or  OD  (points).  The  pull-down  menu  is  used  to  select  the  display 
mode. 
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Viewjnom  allows  the  user  to  select  two  different  windows.  The  first  window  has  the 
data  displayed  as  a  three-dimensional  object,  with  the  vertical  offset  and  color  propor¬ 
tional  to  the  logarithm  of  the  magnitude  of  the  data.  The  second  window  is  similar,  except 
that  it  has  the  color  proportional  to  the  phase  of  the  data.  Each  window  allows  the  surface 
to  be  rotated,  translated,  and  zoomed,  just  as  in  viewjiec.  The  right-hand  mouse  button 
brings  up  the  menu.  The  left  mouse  button  performs  the  transformations.  The  middle 
mouse  button  performs  a  pick  correlation  that  allows  the  user  to  select  a  point  on  the 
matrix.  After  a  matrix  point  has  been  picked  an  information  window  appears  in  the  lower 
left  comer  of  the  console  screen.  This  information  window  describes  the  matrix  and  gives 
information  about  the  selected  point.  Included  in  this  information  are  the  row  and  column 
number  of  the  point.  A  source  code  listing  of  viewjnom  is  in  Appendix  C.  A  sample  MoM 
matrix  display  window  is  shown  in  figure  4. 


Figure  4.  MoM  matrix  2D  display  for  a  four-sided  loop. 

9.0  OBSERVATIONS 

During  this  lED  project  the  authors  learned  a  great  deal  about  visualization  techniques 
and  limitations.  One  observation  is  that  as  more  and  more  features  are  added  to  a  piece 
of  softwEire,  the  user  interface  becomes  the  limiting  factor  in  the  utility  of  the  program. 
Another  observation  has  to  do  with  the  limits  of  hard  copy  devices,  such  as  the  color 
printer  that  was  attached  to  the  SGI.  Although  it  was  a  very  high  quality  printer,  it  was 
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extremely  limited  in  its  ability  to  faithfully  reproduce  the  hue,  saturation,  and  brighmess 
displayed  on  the  console.  Subtle  differences  in  hue  could  not  be  discerned,  and  saturation 
variations  were  washed  out. 

9.1  VIDEO 

In  order  to  properly  document  this  effort,  the  authors  have  proposed  making  a  short 
(10-minute)  video.  This  is  the  only  way  to  demonstrate  the  utility  of  the  codes  Aat  were 
developed.  The  video  is  still  in  production  as  of  the  date  of  this  report. 

9.2  PRODUCTION  CODE 

To  be  used  as  production  codes,  viewjtec  and  viewjnom  need  some  additional  work. 
The  user  interfaces  need  improvement.  NASA’s  Goddard  Space  Flight  Center  has 
developed  a  transportable  applications  environment,  TAE+,  which  provides  a  means  of 
building  a  graphical  user  interface  into  applications.  TAE+  is  designed  to  run  on  the 
Silicon  Graphics  machines. 

There  needs  to  be  a  more  robust  method  of  extracting  NEC  output  data  for  interface 
to  a  NEC  visualization  production  code.  It  has  been  proposed  that  more  CAD-like  capa¬ 
bilities  be  added  to  viewjiec. 

9.3  PC  IMPLEMENTATION 

It  is  not  feasible  right  now  to  port  the  visualization  codes  to  the  PC.  Although  they  are 
written  in  the  C-language,  the  Graphics  Library  is  specific  to  the  Silicon  Graphics 
machine.  However,  SGI  is  attempting  to  make  their  GL  a  graphics  standard  that  will  be 
available  in  the  future  for  PCs.  In  the  meantime,  SGI  is  introducing  a  line  of  desktop 
visualization  machines,  the  Indigos,  with  excellent  graphics  and  CPU  performance  at  an 
exceptionally  low  price,  competitive  with  high-end  PCs. 

9.4  OTHER  CEM  CODES 

Future  work  will  focus  on  bringing  visualization  capabilities  to  other  computational 
codes  used  in  electromagnetic  ship  design.  These  codes  include  the  NEC-Basic  Scattering 
Code  developed  by  Ohio  State  University  and  TSAR,  a  finite  difference  time  domain  code 
developed  by  Lawrence  Livermore  National  Laboratory. 

10.0  CONCLUSIONS 

Several  payoffs  have  resulted,  or  will  result,  from  the  work  done  on  this  lED  project. 

•  Visualization  tools  now  exist  for  performing  analysis  and  validation  at  all  three 
steps  of  a  NEC-MoM  design  procedure. 

•  These  tools  have  greatly  improved  the  utility  of  computational  techniques  used 
in  electromagnetic  analysis  and  are  now  being  used  to  support  topside  ship 
design  projects  for  the  Navy. 
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•  In  combination  with  the  work  done  during  a  related  FY91  lED  project,  the 
visualization  tools  developed  for  this  project  are  supporting  the  efficient  tran¬ 
sitioning  of  computational  electromagnetics  to  high-performance  computing. 

•  High-performance  computing  will  provide  a  huge  increase  in  processing  speed 
and  data  analysis  and  permit  the  more  accurate  modeling  of  ship  topsides. 

•  Since  topside  synthesis  involves  generating  enormous  quantities  of  data,  design 
quality  will  significantly  improve  with  the  enhanced  comprehension  and  compu¬ 
tational  speed  offered  by  visualization  and  high-performance  computing. 

•  This  technology  is  in  the  process  of  being  transitioned  into  the  Electromagnetic 
Engineering  (EME)  project  of  NAVSEA. 

•  Based  on  the  work  done  during  this  lED  project  and  a  related  FY91  lED  pro¬ 
ject,  the  6.2  EMC  project  has  scheduled  projects  in  high-performance 
computing  and  visualization  for  FY94. 
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Appendix  A 

SOURCE  CODE  FOR  filter 


PROGRAM  FILTER 

r 

DIMENSION  XWIREI (3000) , YWIREl (3000) ,ZWIRE1 (3000) , 

SXWIRE2 (3000) , YWIRE2 (3000) , 2WIRE2 (3000) 

DIMENS I(W  XNODEl (3000) , YNODEl (3000) , ZNODEl (3000) , ITAGS(3000) , 
$XNODE2 ( 3000 ) ,YNODE2( 3000) ,2NODE2 (3000), SEGS (3000)  RADS (3000)* 
SSRRAT ( 3000 ) ,RCR( 3000) , PCR( 3000) 

REAL  ICK(3000) ,MCR(3000) 

LOGICAL  CONECT(6000) 

REAL  SEGHIN,SBGHAX,  RAIXIN,  RADMAX,  SRMIN,  SRMAX,  RCRMIN,  RCRMAX, 
SICRMIN, ICRMAX,MCRMIN,MCRMAX.XPCR. RDEL.BRAD, ERAD  LENGTH 
$SLEN,SRAD 

C 

INTEGER  NSEGS, ITAG. IX, ISEG, IRAD, ISR, ICON, IRCR, I  ICR, IMCR 
INTEGER  NWIRES, ISEGS (3000) , JWl RE ( 3000) , JSEG( 3000) 

INTEGER  NNODES, INODE, JNODE, II, 12 
CHARACTER* 12  INFlL 
CHARACTER*T2  LINE 
CHARACTER*!  CH 
C  get  the  fileneme 

WRITE(6,*)  ’Enter  NEC  file  name:  * 

READ(5,10)  INFlL 
10  FORMAT  (A) 

C  open  the  NEC  file 

OPEN(7,  FILE=INFIL) 

C  open  the  temporary  output  file 
OPEN (8,  FILE* 'data . tmp’ ) 

OPEN(9,  FILE-'freq.tmp’) 

C  loop  through  the  NEC  file  reading  the  lines 
20  READ(7, 10, END-50)  LINE 

C  if  its  a  wire  line,  write  it  out 

IF  (LINE(1:2) .EO. ’GW‘)  WRITEt8,*)  LINE(3:) 

IF  (LINE(1:2) .EQ. ’GC*)  WRITEiB,*)  LINE(3:) 

IF  (LINE(1:2) .EO. ’FR*)  WRITE(9,*)  LINEv3:) 

GOTO  20 
50  CLOSE(7) 

CL0SE{8) 

CLOSE (9) 

C  open  up  temporary  frequency  file  and  read  frequency 
OPEN(9,  FILE* * freq, tmp* ) 

FREQ  •  299,8 

READ(9,  *  ,  ENO60)  ISTEP.  ICNT ,  lA,  IB,  FREQ,XSTEP 

60  CLOSB(9) 

C  open  up  temporary  file  and  read  in  wire  lines 
OPEN (7,  FILE»*data . tmp* ) 

C  initialize 

COR-ACOS(0.0)/90.0 
NSEGS  «  0 
NNODES  •  0 

NWIRES  -  0 

C  loop,  reading  data,  also  find  minimum  x,  y,  and  z 
100  READ(7,*, END-200)  ITAG. IX, XI . Y1 , Zl . X2, Y2, 22,X 

NWIRES  «  NWIRES  ♦  I 
ITAGS (NWIRES)  =  ITAG 
ISEGS (NWIRES)  -  IX 
XWIREI (NWIRES)  -  XI 
XWIRE2 (NWIRES)  =  X2 
YWIREl  (NWIRES)  -  Y1 
YWIRE2  (NWIRES)  -  Y2 
ZWIREl (NWIRES)  »  Zl 
ZWIRE2 (NWIRES)  -  Z2 
IF  (NSEGS. EQ.O)  THEN 


XHIN  - 

XI 

YMIN  » 

Yl 

ZMIN  - 

Zl 

SEGMIN 

•  10000.0 

RADHIN 

-  10000.0 

SRMIN  > 

=  10000.0 

RCRMIN 

-  10000.0 

ICRMIN 

-  10000,0 

MCRMIN 

•  10000.0 

XMAX  - 

XI 

YMAX  - 

Yl 

ZMAX  - 

Zl 

SEOtAX 

-  0.0 

RADNAX 

-  0.0 

SRMAX  » 

'  0.0 

RCRMAX 

-  -10000.0 

ICRMAX 

-  -10000.0 

MCRMAX 

-  0.0 

END  IF 

IF 

(Xl.LT.XMIN) 

XMIN-Xl 

IF 

(Yl.LT.YMIN) 

YMIN-Yl 

IF 

(Zl.LT.ZMIN) 

ZMIN-Zl 

IF 

(XI .GT.XMAX) 

XMAX-Xl 

IF 

(Y1.GT.YMAX) 

YMAX-Yl 

IF 

(Zl.GT.ZMAX) 

ZMAX-Zl 

IF 

<X2.LT.XMIN) 

XMIN-X2 

IF 

(Y2.LT.YMIN) 

ymin-y: 

IF 

(Z2.LT.2MIN) 

ZMIN-22 
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IF  (X2.GT.XMAX)  XMAX«X2 
IF  (Y2.GT.VMAX)  YMAX-y2 
IF  (Z2.GT.ZMAX)  ZMAX-Z2 
IF  (X.GB. 0.000001)  THEN 
DO  120  I-l,IX 
NSEGS  -  NSEGS+1 
NNODES  -  NNODES+2 
JWIRE(NSEGS)  -  NWIRES 
JSEG(NSEGS)  -  I 

XNODBl(NSEGS)  -  XI  ♦  { I-l )  MX2-X1 )  /  IX 
YNODBl (NSEGS)  •  Y1  ♦  { I- 1 ) * (Y2- Yl ) / IX 
ZNODBl (NSEGS)  »  Z1  ♦  ( I-l ) • (Z2-Z1 ) / IX 
XNODE2(NSEGS)  -  XI  ♦  IMX2-X1)/IX 
YNOOE2(NSEGS)  -  Y1  ♦  I*(Y2-Y1)/IX 
ZNODE2(NSEGS)  -  Zl  ♦  IMZ2-Z1)/IX 

SEGS(NSEGS)  -  SQRT ( (Xl -X2) *  *  2+ (YI-Y2) **2+ (21-Z2) • •2)/IX 
RADS (NSEGS)  -  X 

SRRAT(NSEGS)  «  SEGS (NSEGS) /RADS (NSEGS) 

CONBCT (NNODES)  »  .FALSE. 

CONECT(NNODES-l)  -  .FALSE. 

IF  (SEGS (NSEGS) .LT.SBGMIN)  SEGMIN-SEGS (NSEGS) 

IF  (RADS (NSEGS) .LT.RADMIN)  RADMIN-RADS (NSEGS) 

IF  (SRRAT(NSBGS) .LT.SRMIN)  SRMIN-SRRAT (NSEGS) 

IF  (SEGS(NSEGS) .GT.SEGMAX)  SEGMAX-SEGS (NSEGS) 

IF  (RADS (NSEGS) .GT.RADMAX)  RADMAX-RADS (NSEGS) 

IF  (SRRAT(NSEGS)  .GT.SRMAX)  SRMAX**SRRAT (NSEGS) 

120  CONTINUE 
ELSE 

READ (7,*, END-200)  11 , 12. RDEL, BRAD, ERAD 
DO  140  I-l, IX 
NSEGS  -  NSEGS* 1 
NNODES  «  NNODES*2 
JWI RE (NSEGS)  -  NWIRES 
JSBG(NSEGS)  -  I 
IF  (I.GT.l)  THEN 

XI  •  XNODE2(NSEGS-l) 

Yl  -  YNODE2 (NSEGS- 1) 

Zl  -  ZNODE2 (NSEGS- 1) 

END  IF 

LENGTH  «  SQRT((X1-X2)**2+(Y1-Y2)**2*(Z1-Z2)**2) 

IF  (I.EO.l)  THEN 

IF  (ABS(l-RDEL) .GT.O.OOOl)  THEN 
SLEN  •  LENGTHMl-RDBL)/(l“RDEL**rX) 

ELSE 

SLEN  •  LENGTH/ IX 
END  IF 

ELSE  IF  (I.EO.IX)  THEN 
SLEN  •  LENGTH 
ELSE 

SLEN  •  ROEL*SLEN 
END  IF 

SRAD  -  BRAD  +  ( I-l ) * (ERAD-BRAD) / ( IX-1 ) 

XNQDEI (NSEGS)  «  XI 
YNODEl (NSEGS)  -  Yl 
ZNODEl (NSEGS)  -  Zl 

XNODE2(NSEGS)  -  X1+ (SLEN/LENGTH)  MX2-XI) 

YHODB2 (NSEGS)  «  Yl ♦ (SLEN/LENGTH) ♦ (Y2-YI ) 

ZNOOE2 (NSEGS)  -  21* (SLEN/LENGTH) * (Z2-Z1 ) 

SEGS (NSEGS)  -  SLEN 
RADS (NSEGS)  «  SRAD 
SRRAT (NSEGS)  «  SLBN/SRAD 
CONBCT (NNODES)  -  .FALSE. 

CONECT(NNODES-l)  -  .FALSE. 

IF  (SEGS (NSEGS) .LT.SEGMIN)  SEGMIN-SEGS (NSEGS) 

IF  (RADS(N.SEGS)  .LT.RADMIN)  RADMIN-RADS  (NSEGS) 

IF  (SRRAT(NSBGS) .LT.SRMIN)  SRMIN-SRRAT (NSEGS) 

IF  (SEGS (NSEGS) .GT.SEGMAX)  SEGMAX-SEGS (NSEGS) 

IF  (RADS (NSEGS) .GT.RADMAX)  RADMAX-RADS (NSEGS) 

IF  (SRRAT(NSEGS) .GT.SRMAX)  SRMAX-SRRAT (NSEGS) 

140  CONTINUE 
END  IF 
GOTO  100 
200  CLOSE  (7J 

C  se«  which  segments  ere  connected 
DO  300  I-l, NSEGS 
INODE  -  (I*2)-l 

IF  (ZNODEl (I) .LE.RADS(I) )  CONBCT(INODE)-.TRUE. 

IF  (2NODE2(I).LE.RADS(I))  CONE(rr(INODE*l)-.TRUB. 

DO  350  J-1, NSEGS 
IF  (J.NB.I)  THEN 
JNODE-(J*2)*-l 

IF  (.NOT.CONECT(INODE))  THEN 
DIST  -  SQRT((ZNODEl(J)-ZNODEl(I))**2* 

S  (YNODEl (J) -YNODEl (I) ) ••2* (XNODEl ( J> -XNODBl ( I ) ) *  *  2) 

IF  (DIST.LE.RADSd)  )  THEN 
CONECT(INODE)  -  .TRUE. 

CONB(rr(JNODE)  -  .TRUE. 

END  IF 

DIST  -  SQRT((2N0DE2(J)-ZN0DE1(I))**2* 

S  (YN0DE2(J) -YNODBl (I) )**2*(XMODE2{J)-XNODEl (I))**2) 

IF  (DIST.LE.RADS(I))  THEN 
CONECT(INODB)  -  .TRUE. 

CONECT(JNODE*l)  •  .TRUE. 
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END  IF 
END  IF 

IF  ( .NOT.CONECT(INODE+l) )  THEN 

DIST  «  SQRT((2NODEl(J)-ZNODE2(I))**2*' 

$  (YNODEl (J) -YNODE2(I) )**2+ (XNODEl ( J) -XNOOE2 ( 1) ) * *2 ) 

IF  (DIST.LE.RADS(I))  THEN 
CONBCT( INODE* 1)  -  .TRUE. 

CONECT(JNODE)  -  .TRUE. 

END  IF 

DIST  -  SQRT((ZNODE2(J)-ZNODE2(I))**2* 

$  ( YNODE2 ( J ) «  yNODE2 ( I ) ) *  *  2  * ( XNODF2 { J ) -XNODE2 { I ) ) *  *  2 ) 

IF  (DIST.LE.RADSd))  THEN 
CONECT(INODE*I)  »  .TRUE. 

CONECT{JNODE+1 )  «  .TRUE. 

END  IF 
END  IF 

IF  <CONECT(INODE) .AND,CONECT(INODE+l))  GOTO  300 
END  IF 
3S0  CONTINUE 
300  CONTINUE 
C  calculate  shifts 

XSHIFT  -  {XMAX  ♦  XMIN)/2.0 
YSHIFT  *  (YMAX  +  YHIN)/2.0 
ZSHIFT  -  (ZMIN) 

C  calculate  scale  factor 
XDIF  »  XMAX-XMIN 
YDIF  -  YMAX-YMIN 
ZDIF  •  ZKAX-ZMIN 
DIFMAX  «  XDIF 

IF  (YDIF. GT. DIFMAX)  DIFMAX=YDIF 
IF  (ZDIF. GT. DIFMAX)  DIFMAX«ZDIF 
SCALE  =  4.0/DIFMAX 
C  get  currents  if  necessary 
C  see  if  there  is  an  output  currents  file 

WRITE(6,*)  ’Is  there  an  output  currents  file?  (Y:Y€3,  NtNo) : * 
READ(5,10)  CH 

IF  (CH.EQ. ’N'.OR.CH.EQ. ’n')  THEN 
DO  320  I*1,NSEGS 
RCR(I)  -  0.0 
ICR(I)  -  0.0 
MCH(I)  «  0.0 
PCR(I)  »  0.0 
320  CONTINUE 

RCRMIN  »  0.0 
RCRMAX  -  0.0 
ICRMIN  -  0.0 
ICBMAX  *0.0 
MCRMIN  -0.0 
MCRMAX  »  0.0 
ELSE 

C  get  the  filename 

WRITE(6,*)  'Enter  output  currents  filename:  ' 

READ(S,10)  INFIL 
C  open  the  currents  file 

OPEN(T,FlLE«INFIL) 

C  loop  through  reading  currents  (magnitude,  phase) 

RBAD(7,10)  line 
READ(7,10)  LINE 
DO  340  I-1,NSEGS 

READ(7,M  MCR(n,PCR(I) 

RCR(I)  -  MCR(I)^COS(CDR*PCR(I)) 

ICRd)  -  MCR(I)*SIN(CDR*PCR(I)> 

IF  (MCR(I) .LT. MCRMIN)  MCRMIN=MCR(I) 

IF  (MCR(I) .GT. MCRMAX)  MCRMAX*MCR(I> 

IF  (RCR(I) .LT. RCRMIN)  RCRMIN=RCR(I) 

IF  (RCR{I) .GT. RCRMAX)  RCRMAX=RCR ( I ) 

IF  (ICR(I) .LT. ICRMIN)  ICRMIN»ICR ( I ) 

IF  (ICR(I).GT.ICRMAX)  ICRMAX»ICR( I ) 

340  CONTINUE 
CL0SE(7) 

END  IF 

C  new  print  it  all  out 
C  get  the  filename 

WRITECS,*)  'Enter  output  file  name:  ' 

READ(5,10)  INFIL 
C  open  the  output  file 

OPEN(8,  FILE-INFIL) 

WRITB(8,*)  FREO 
WRITE(8,*)  NWIRES,  NSEGS 
WRITE (8, M  XMIN,  XMAX 
WRITE(8,*)  YKIN,  YMAX 
WRITE{8,*)  ZMIN,  ZMAX 
WRITE(8,*)  XSHIFT,  YSHIFT 
WRITB(8,*)  ZSHIFT,  SCALE 
WRITE(8,*)  SEGMIN,  SEGMAX 
WRITEO,*)  RADMIN,  RADMAX 
WRITE(8,*)  SRMIN,  SRMAX 
WRITE(8,*)  RCRMIN,  RCRMAX 
WRITE(8,M  ICRMIN,  ICRMAX 
WRITE(8,*)  MCRMIN,  MCRMAX 
DO  360  I-l, NWIRES 

WRITE  (8, 902)  ITAGS(I).ISEGS(I),XWIREl(n,YWIREUl),ZWIREUI). 
SXWIRE2(I) .YWIRE2(I) ,ZWIRE2(I) 


A-3 


360  COWTINUE 

902  FORMATC  ',16,14,6210.4) 

DO  400  I-1,NSEGS 
INODE-(I*2)-l 

IF  (SEGMAX.EQ.SEGMIN)  THEN 
ISEG  -  0 
ELSE 

tSEG  -  (SEGS(Ii-SEGMrN)/(SEGHAX-SEGMIN)*7 
IF  (SEGS(I) .EQ.SEGMAX)  ISEG  -  6 
END  IF 

IF  (RADMAX.Ee.RADMIN)  THEN 
IRAD  -  0 
ELSE 

IRAD  -  IRADS(l) -RADMIN)/ (RADHAX-RADMIN)*? 

IF  (RADS(I)  .EQ.RADMAX)  IRAD  -  6 
END  IF 

IF  (SRMAX.EQ.SRMIN)  THEN 
ISR  -  0 

ISR  -  (SRRAT(I)-SR14IN)/(SRMAX-SRMIN)*7 

IF  (SRRAT(I) .EQ.SRMAX)  ISR  -  6 
END  IF 

IF  (RCRMAX.EQ.RCRMIN)  THEN 
IRCR  -  0 
ELSE 

IRCR  -  (RCR(I) -RCR.NIN)/ (RCRMAX-RCRMIN)*T 
IF  (RCRII) .Eg.RCRMAX)  IRCR  -  6 
END  IF 

IF  (ICRMAX.EQ. ICRMIN)  THEN 
IICR  -  0 
ELSE 

IICR  -  (ICR(I)-ICRmN)/(ICRMAX-ICRMIN)‘7 
IF  (ICR(I) .EO. ICRMAX)  IICR  =  6 
END  IF 

IF  (MCRKAX.EQ.MCRMIN)  THEN 
IMCR  -  0 
ELSE 

IMCR  -  (MCR(I) -MCRMIN)/ (MCRMAX-MCRMIN)*7 
IF  (MCR(I) .EQ.MCRHAX)  IMCR  -  6 
END  IF 

C  put  XPCR  between  0  to  360 
XFCR  -  PCR(I) 

410  IF  (XPCR.LT.O)  THEN 

XPCR-XPCR«360.0 
GOTO  410 
END  IF 

411  IF  (XPCR. GE. 360.0)  THEN 

XPCR-XPCR-360.0 
GOTO  411 
END  IF 

IF  (CONECT(INODB)  .AND.CONECTlINODE+U)  ICON  =  0 
IF  (CONECT(INODE).AND.  (.NOT.CONECTdNODE*!)))  ICON  =  2 
IF  (  (.NOT.CONECT(INODE)  )  .AND.CONECT(INODEH) )  ICON  =  2 
IF  ((.NOT.CONECT(INODE)).AND.(.NOT.CONECT(INODE*l)))  ICON  »  4 
WRITE{9,903)  JSEG(I),JWIRE(I),ISEG. IRAD, ISR, 

$  ICON, IRCR, IICR, IMCR, XPCR 

903  FORMATC  ’,15,16,713,210.4) 

WRITE (8, 905)  SEGS(I) .RADS(I) ,SRRAT(I) ,RCR(I),ICR(I),MCR(I) 

'6  FORMATC  ’,6210.4) 

WRITE(e,904)  (XNODEKD),  (YNODEKD), 

S(2N0DE1 (I) ) , (XN0DE2(I)) , 

S(VNODE2(I)), (ZNODE2(I)) 

904  FORMATC  ’,6210  4) 

400  CONTINUE 

CLOSE(8) 

END 
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Appendix  B 

SOURCE  CODE  FOR  viewjtec 

View  n«c .c 

vt9w  rcc  displays  t  n«c  wit«  Th«  pol«cvi«w(| 

trtnsf otMst  i<4)  is  chan)]«d  with  aovis*  input  wh«>n 
i0ft  hutrf>n  19  HbKI  down.  Wire  p«r«a«teta  are 

usin^  fha  aiidll#  saih#  Menn9  are 

uf  j  the  figh?  aosjsa 

i  y*  4  j^‘-  1^'H.  SJlicrn  Giajhira,  Inc,  Ri^jhts  tteaerved. 

‘  rciyri-jh?  f nt  rifjhta  and  Habiliry  rtp.at i-  n. 

r.  t:  Tahnital  R.1u-a?iv.s  ^‘uijina  D«»veli  {^er^ 


•  ji  \j,  h» 

iU'-  di  la  '  Ti  '  !#*  i  •  .  h  ' 
I  if  .  ^  t#  •  9*  J  l  ’  ?.  • 

•  .-.  ,  ;  !a  *  h,  K 


!  »  .  'J%  f  t '  I  »  yjw» 


,  ,1  ini?  iai  kza  •  ; 

>.  i1  ■air.ivM’15  ; 

V  i1  prin?h4?9{{  r.;}  itiirti; 

V  k-l  fiav 'rgb!  f  »  '<a?  hs.* 

V  .-1  lraw3'car»e’Shrtt  'Jut*  shoi?  dx? .  -set?  dxz, 

dyr,  shrf»  -iy?,  ?hnr?  dy:, 
ir»r.d  data  chci:»,  i  ni  ? rana^ch^ica,  lonq  ibllnk« 

t'uf?ar!l,  i-r.:i  hi*  a,  B«>^laan  picRin-j,  lon^  do^win)? 

'•it  aaea  .void ' . 

V  i 3  initethiv -  1 1) i 
V.'. 4  3  ini?efh(v-.>j1).‘ 
veil  inltfiii  v*  id} ; 

y  4  3  init»nd#I  tv:;l'3> ; 

veil  «>aaodaI  Uc-n  ]  data  rhoice.  iorij  iblinki; 


flLC  •|eofil«; 
rtU  •atvUle; 

Flic  *etpfne; 
f:M  «aF«fila; 

Fits  *eppflle; 

FILE  •FldfiUj 
f'tt  file,naM{^y)|; 
tar 

r.ac  ettt^fllei.'^I; 

■tar  atp^flleU'^J? 
hat  apai^f  ilaf-SJ ; 
hat  epf*  fllatrS); 

'bar  rid_rilei:sw 
9*ati:  char  deoH  *  ".law*; 
i*a?lr  thar  eteil  *  •.«?«•; 
i»a*.ic  ‘hat  e^rll  •  “.atp"; 

9?a*if  'har  eF**!!  •  •.efta'*; 
ftati^  i:hat  eprJl  •  •.epp"; 
ata?3c  char  fllH  *  •.fid"; 
ir?  sthreth  •  0; 

in?  frhreah  • 

if*  nwirea,  naaga.  »aa-j«  I .  ?a«  j  I  ,  iwirellOOOl,  data  (3000]  f7|  ; 

irt?  l?a<39(  33001 ; 

9*aTic  Boolean  iraw^axea  •  TRUE; 

fi:*a?  b^.r<^vdeI3'300iT3I,  eivin<r)deU'^''^3}  (  M  .  b^jnwire  (  iOOOM  3| »  endwlreOOOOJ  1 3| ; 
:ice*  lutailn,  AxiaaR,  ayvin.  ayiaaM,  aznir.,  aissaa.  sacale.  Ndata(3000]; 
floM*  xshift»  yahlft,  ishlft,  freq,  harp',  hero; 

flat  aeganoOOl  «rads|  3C001 ,9rrat  |  3000) ,  rcr  1 30001 ,  \cr  (30001  ^mer  [3000] ; 
flcMit  aeqetin.  aeqaea,  radaiin,  ra<baaK.  arBin,  sraax,  rermin,  rcraax* 

ircBio,  icratax.  »craiin,  acraiak,  ethaiax.  ephsax.  fldmax,  zfldMx; 

if#-'!  rqfecolfili 


static  fl  iat  tx.  *f.  ts; 

/•  wifvirM  id's 

L  aeq  Win,  ra-i^wm,  rat^win,  ,  real^win,  laaq^win. 

•aq_,win,  ph  win,  eth  win,  ep*h^win,  dale_^win,  cur  win, 
fld^win,  sfld^win; 

static  char  teq^win^naaiel  I  •  "Sed^efitat  iryn"; 
a*atic  char  rad_win  naatelj  •  "Wife  radius**; 
static  char  rat  win  na»e(]  ’•Seqiaen?  to  tadiua**; 

•tattc  char  «‘<>n  win^naAeU  •  •r/winect  lona"; 
sta?i-  char  real^wln  nasMall  “  "Real  cottent**; 
static  <-haf  laaf^wln  na»e  I  ]  •  "laaqinafy  .-uriefst**; 
s*a*«’  that  Mq  win^na**!)  *  •Majnitode  r,t  current  ■*; 
s*a»ir  rhaf  ph  win  naSMti  •  •Fhase  i  f  iirren?**; 

•*atic  char  ety.^wln  naaeU  *  **S  Iheta^i 

at*t4c  char  ef-t  win  na»e|  |  -  '•R  phi**; 

atatir  char  da? a^win_na»e| )  •  *'Wire  Raraseteis**; 

statit  chat  f !d_win  naMi }  •  "Wear  field**; 

static  char  rf Id  wtn_na»#f I  »  •2*cf«p  of  Near  field**; 

l’4ii  loaax.  ysMM,  zwal; 


Reproduced  From 
Best  Available  Copy 


B-l 


/*  initialize 

*  The  initialize  subroutine  positions  the  windc»<  and  specifies 

*  its  future  constraints.  Graphics  configurations  are  set  and 

*  the  event  queue  is  initialized. 

♦/ 

void 

initialize (void) 

( 

printf(”  Fnter  tile  name;  "); 
scanf("ls",  file^name) ; 
strcpy (geo_file, £ile_name) ; 
strcpy (etm_file, file_name) ; 
strcpy (etp_file, file^name) ; 
strcpy  (epin_file,  file_name) ; 
strcpy (epp_f ile, file_narae) ; 
strcpy ( f ld_f ile, file_name) ; 
street (geo_file, geo) ; 
street (etm_file, etm) : 
street (etp_file, etp) ; 
street  (ep(n_file,epa) ; 
street (epp_file,epp) ; 
streat (fld_file, fid) ; 

if((9eofile  *  fopen (geo_file, "r") )  «■  (FILE*)0) 

I 

printf ("XnXtXt  File  Is  could  not  be  opened,  \n**,  geo^file); 
exitO  ; 

1 

xmax  »  getgdesc (GD_XPHAX) ; 
ynax  »  getgdesc  (GD_^YPMAX) ; 
zval  «  getgdesc (GD_2MAX) ; 

minsize (xaax/ 10,  ymax/lO); 
maxsize (xfflax-200,  yniax-160); 
keepaspect (xmax,  ymax); 
seg__win  »  winopen  ('*") ; 
wintitle(seg_win_name) ; 

winconstraints ( ) ; 
winset (seg^win) ; 
zbuffer(TRUE); 
doublebuffer  0; 

RGBmode  ( ) ; 
gconfigO  ; 

shademodel (FLAT) ; 

Raaode(KVIEWING) ; 

qdevice(LEFTMOUSE); 
qdevlce(MIDDLEMOUSE) t 
qdevice(RIGHTMOt;SE) ; 
qdevice(ESCKEY); 
tie(LEFTMOUSE,  MOUSEX,  MOUSEY); 


I  /♦  end  initializeO  */ 


/*  main 

*  main  calls  initialize,  then  goes  into  a  loop.  Within  the  loop 

*  drawscene  is  called  and  events  in  the  event  queue  are  processed 

*  until  user  exits  (by  pressing  ESCAPE  or  through  the  window  manager) . 

*  Mouse  input  is  passed  to  drawscene  when  the  left  mouse  txjtton 

*  is  held  down. 

*/ 

void 

main (void) 

( 

Boolean  exitflag  >  FALSE; 

Boolean  redraw^needed  »  TRUE; 

Boolean  draw_seg  ■  TRUE; 

Boolean  draw_rad  *  FALSE; 

Boolean  draw_rat  •  FALSE; 

Boolean  draw_con  ■  FALSE; 

Boolean  draw^real  •  FALSE; 

Boolean  draw_imag  *  FAI^E; 

Boolean  draw_mag  •  FALSE; 

Boolean  draw_ph  ■  FALSE; 

Boolean  draw_eth  ■  FALSE; 

Boolean  draw_eph  *  FAI«SE; 

Boolean  draw_data  ■  FALSE; 

Boolean  draw_^fld  ■  FALSE; 

Boolean  draw_zfld  »  FALSE; 

short  value; 
short  buffer (BUFSIZE) ; 
long  hits  •  0; 

Boolean  picking  «  FALSE; 
short  deltax,  deltay; 

long  dev; 
long  choice; 

long  data_choice,  trans^choice,  iblink; 
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long  m«nu,  datamenu,  transmenu,  geomenu,  curmenu, 
nfieldmenu,  ffieldmenu; 

short  X,  oldx;  /*  mouse  movement  */ 

short  y,  oldy; 
short  dxr,  dxt,  dxz  -  0; 
short  dyr,  dyt,  dyz  =  0; 
static  long  whitecoid  *  (  255,  255,  255  J; 
static  long  bluecol(]  *  i  0,  0,  255  ); 
static  long  dbluecoK]  «  I  0,  0,  128  ); 

initialize ( ) ; 
data_choice  *1; 
trans  choice  «  1; 
iblink  »  0; 
tx  «  0; 
ty  »  0; 
tz  -  0; 

/*  make  the  geometry  menu  */ 

geomenu  *  defpup('*Choose  Geometry  Data  Window  ttl** 

"Segmentation  txllWire  radius  lx2ISegment  to  radius  ratio  tx3l** 

"Wire  connections  tx4"); 

/*  make  the  currents  menu  •/ 

curmenu  *  defpup( "Choose  Currents  Data  Window  It  I" 

"Real  ctxnponent  of  current  1x5 1" 

"Imaginary  component  of  current  1x6 (Magnitude  of  current  1x71" 

"Phase  of  current  1x8"} ; 

/*  make  the  far  fields  menu  ♦/ 

ffieldmenu  »  defpup( "Choose  Far  Fields  Data  Window  Itl" 

"E-theta  component  1x9 1  E-phi  component  1x10"); 

/*  make  the  near  fields  menu  */ 

nfieldmenu  -  defpupC’Choose  Near  Fields  Data  Window  Iti" 

"2-component  of  E-normal  1x11 (Total  E-notmal  1x12"); 

/♦  make  the  data  menu  •/ 

datamenu  ■  defpup{"Choose  Data  Window  to  Open/Close  %t|" 

"Geometry  ImICurrents  lm(Near  Fields  ImiFar  Fields  Im", 
geomenu,  curmenu,  nfieldmenu,  ffieldmenu); 

/♦  make  the  transformation  menu  ♦/ 
tranamenu  »  defpup("Choose  Transformation  Itl" 

"Rotate  lxl3lTranslate  IxllIZoom  IxlSIReturn  to  initial  state  1x16"); 

/♦  make  the  main  menu  ♦/ 

menu  ■  defpup("Menu  It (Data  Window  ImITransformation  1ml" 

"Close  Current  Window  IxlT (Toggle  Axes  On/Off  1x18", datamenu, transmenu) 

while  (exitflag  »  FALSE) 

if  (redraw  needed  »»  TRUE) 

I 

if  (draw_seg) 

drawscene(dxr,  dxt,  dxz,  dyr,  dyt,  dyz, 

1,  tran3_choice,  iblink, 
buffer,  hits,  picking,  3eg_win) ; 

if  (draw_^iad) 

drawscene (dxr ,  dxt,  dxz,  dyr,  dyt,  dyz, 

2,  trans^choice,  iblink, 
buffer,  hits,  picking,  rad_win) ; 

if  (draw_rat) 

drawscene(dxr,  dxt,  dxz,  dyr,  dyt,  dyz, 

3,  ttans_choice,  iblink, 
buffer,  hits,  picking,  rat_win); 

if  {draw_con) 

drawscene(dxr,  dxt,  dxz,  dyr,  dyt,  dyz, 

4,  trans_choice,  iblink, 
buffer,  hits,  picking,  con_win); 

if  (draw_real) 

drawscene(dxr,  dxt,  dxz,  dyr,  dyt,  dyz, 

5,  trans_choice,  iblink, 
buffer,  hits,  picking,  ceal_win); 

if  (draw_imag) 

drawscene (dxr ,  dxt,  dxz,  dyr,  dyt,  dyz, 

6,  trans_choice,  iblink, 
buffer,  hits,  picking,  imag_win); 

if  (draw_mag) 

drawscene(dxr,  dxt,  dxz,  dyr,  dyt,  dyz, 
f ,  trans_choice,  iblink, 
buffer,  hits,  picking,  mag_win) ; 

if  (draw_ph) 

drawscene (dxr ,  dxt,  dxz,  dyr,  dyt,  dyz, 

8,  trans^choice,  iblink, 
buffer,  hits,  picking,  ph  win); 

if  (draw_eth) 

draw3cene(dxr,  dxt,  dxz,  dyr,  dyt,  dyz, 

9,  trans_choice,  iblink, 
buffer,  hits,  picking,  ethwin); 

if  (draw_eph) 

drawscene (dxr ,  dxt,  dxz,  dyr,  dyt,  dyz, 

10,  tran3_choice,  iblink, 
buffer,  hits,  picking,  eph_win); 

if  (draw_2fld) 

drawscene (dxt ,  dxt,  dxz,  dyr,  dyt,  dyz, 

11,  trans_choice,  iblink, 
buffer,  hits,  picking,  zfld_win); 

if  (drawfld) 


drawscen* (dxr,  dxt,  dxz,  dyr,  dyt.  dyz, 

12,  trans_choice,  iblink, 
buff«r,  hits,  picking,  fldjwin); 

if  (iblink  •-  0)  redraw^needed  »  FALSE; 
iblink  •  -iblink; 

1 

while  (  (exitflag  »  FALSE)  <c 

(qtestO  II  (redraw  needed  FALSE)}  ) 

{ 

dev  ■  gread  (cvalue); 
if  (dev  «  ESCKBY) 

I 

if  (value  -»  0) 

exitflag  *  TRUE; 

) 

else  if  (dev  —  REDRAW) 

( 

reshapeviev^it  ( ) ; 
redraw_needed  *  TRUE; 

1 

else  if  (dev  INPUTCHANGE) 

I 

cur  win  »  value; 

I 

else  if  (dev  RIGHTM<XISE  44  value  1) 

I 

choice  >  dopup(menu); 
if  (choice  >-  1) 

( 

if  (choice  ■■  1) 

( 

draw^seg  •  fdtaw^seg; 
if  (draw  seg) 

( 

minsize (xnax/ 10, ymax/ 10) ; 
maxsi2e(xnax-200,ymax-160) ; 
keepaspecc (xmax, ymax) t 
seg^win  ■  winopen  (*"*) ; 
wintitle (seg_win_naae) ; 
winconstraintsO ; 
winset(seg  win); 

2buf/er(TRU£)/ 

doublebufferO; 

RGBmode ( ) ; 
gconfigO; 
shademodel ( FLAT) ; 
mffiode{KVI EWING); 

) 

else  winclose(seg  win); 

I 

if  (choice  -*  2) 

( 

draw_rad  -  'draw_rad; 
if  (draw_rad) 

( 

minsize (xmax/ 10, ymax/ 10) ; 
maxsize (xmax-200, ymax- 160) ; 
keepaspect (xmax, ymax) ; 
rad_win  «  winopen ("") ; 
wintitle (rad_win_name) ; 
winconstraints ( ) ; 
winset (rad^win) ; 
zbuffer (TRUE) ; 
doublebuffer ( ) ; 

RGBmode ( ) ; 
gconf ig ( } ; 
shademodel(FLAT) ; 
mBode(HVI  EWING); 

1 

else  winclose(rad_^win) ; 

I 

if  (choice  --  3) 

( 

draw^rat  ■  !draw_rat; 
if  (draw_rat) 

( 

minsize (xmax/ 10, ymax/ 10) ; 
maxsize (xmax-200, ymax- 160) ; 
keepaspect (xmax, ymax) ; 
rat_win  ■  winopen  (*•") ; 
wintitle(rat_win_naree) ; 
winconstraints () ; 
winset (rat_win) ; 
zbuffer (TRUE); 
doublebuffer () ; 

RGBmode ( ) ; 
gconf ig ( ) ; 
shademodel (FLAT); 
nmode(KVIEWlNG); 

I 
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else  winclose trat_win) ; 


if  (choice  -«  4) 

i 

<ltaw_con  »  !draw_con; 
if  (draw  con) 

I 

minsize  (xmax/ 10,  yinAK/lO)  ; 
maxsize (Kmax-200, ymax-160) ; 
keepaspect (xmax, ymax) ; 
con  win  “  winopen  (**") ; 
wintitle (con^winname) ; 
winconstraintsO  ; 
winset (con_win) ; 
zbuffer(TRaE); 
doublebuffer (} ; 

RGBmode ( ) ; 
gconf ig {) ; 
shademodel (FLAT) ; 
ininode(MVI  EWING): 

} 

else  winclose  (c'^n_win> ; 

1 

if  (choice  »=  f>) 

I 

draw_real  *  !dtaw_real; 
if  (draw_real) 

( 

tninsize  (xmax/ 10,  ymax/ 10) ; 
maxsize (xmax-200, ymax- 160) ; 
keepaspect (xmax, ymax) ; 
real^win  *  winopen<****) ; 
wintitle (real_win_name) ; 
winconstraints ( ) ; 
winset { real_win) ; 
zbufferlTRUE): 
doublebuf fee () ; 

RGBmode ( ) : 
gconf Ig ( ) ; 
shademodel  (Fl-AT)  ; 
mBode(HVIEWING) ; 

1 

else  winclose(real_win) ; 

I 

if  (choice  6) 

I 

draw  imag  *  !draw_imag; 
if  (draw  imag) 

( 

minsize (xmax/ 10, ymax/ 10) ; 
maxsize (xmax-200, ymax- 160) ; 
keepaspect (xmax, ymax) ; 
imag_win  »  winopen  (*•••) ; 
wint itle(imag_win_name) ; 
winconstraints ( ) ; 
winset (imag  win); 

Zb'if fee  (TRUE); 
doublebuf  fer ( ) ; 

RGBmcde ( > / 
gconf ig( ) ; 
shademodel (FLAT) ; 
mmrde{MVI EWING) ; 

I 

else  winclose (imag_win) ; 

I 

if  (choice  *-  T) 

( 

draw  mag  -  'djaw_mag; 
if  (drawmag) 

I 

minsize (xmax/ 10, ymax/ 10) ; 
maxsize (xmax-200. ymax- 160) ; 
keepaspect (xmax, ymax) ; 
mag  win  *  winopenC***); 
wintitle (magwinname) ; 
winconstraints () ; 
winset  (mag__win) ; 
zbuffer(TRUE); 
doublebuffer ( ) ; 

RGBmoJe ( ) ; 
gconfigO  ; 
shademodel (FLAT) ; 
mmole (MVI EWING); 

I 

else  wine  lose  {m^*g  win); 

) 

if  (choice  -■=  8) 

\ 

draw_ph  *  Idrawph; 
if  (drawph) 

( 

minsize (xmax/ 1 0, ymax/ 10) ; 
maxsize (xmax -200, ymax -160); 
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k««p«sp«ct (xmax, y«ax) ; 
ph_win  ■  winop«n<****) ; 
wintitla (ph_win_name) ; 
winconstraints ( ) ; 
winset (ph_win) ; 
zbuf far (TRUE); 
doublabuffer () ; 

RGBmode ( ) ; 
gconfigO  ; 
shadamodeKFLAT); 

Bmode(MVI  EWING); 

) 

elsa  winclo3a(ph_win) ; 

) 

if  (choica  -■  9) 

{ 

draw_ath  •  !draw^ath; 
if  (draw^ath) 

{ 

minaiza (xmax/ 10,ymax/10) ; 
maxaiza (xmax-200, ymax>160) ; 
kaapaapact (xmax, ymax) ; 
ath^win  •  winopen(“**) ; 
wintitla(ath_win_nama) ; 
winconatrainta ( ) ; 
winaet (eth_win) ; 
zbuffar(TRU£) ; 
doublabuffer ( ) ; 

RGBmoda ( } ; 
gconfigO  ; 

ahademodel (GOURAUO) ; 
mmoda(HVI  EWING); 

» 

elae  wincloae(eth  win); 

) 

if  (choice  10) 

I 

draw^eph  ■  !draw_^eph; 
if  (draw^eph) 

I 

minaiza (xmax/ 10, ymax/ 1 0) ; 
maxaiza (xmax-200, ymax- 160) ; 
kaapaapact (xmax, ymax) ; 
aph^win  ■  winopani"") ; 
wintitia(aph^win_nama) ; 
winconatrainta (); 
winaet (aph  win); 
zbuf far (TRUE); 
doublabuffer () ; 

RGBmoda ( ) ; 
gconfigO ; 

ahademodel (GOURAUD) ; 
aaBoda(HVl  EWING) ; 

I 

alaa  wincloaa(aph_win> ; 

» 

if  (choica  ■«  II) 

I 

draw_zfld  -  'draw_zfld; 
if  (draw_zfld) 

( 

minaiza (xmax/ 10, ymax/ 10) ; 
maxaiza (xmax-200, ymax- 160) ; 
kaapaapact (xmax, ymax) ; 
zfld_win  •  winopan {•••*)  ; 
wintitla(zfld_win_nama) ; 
winconatraintaO ; 
winaat (zfld^win) ; 
zbuffar(TRUE) ; 
doublabuffer  O ; 

RGBmoda ( ) ; 
gconfigO  ; 
ahademodel (FLAT) ; 
mmoda(MVI EWING); 

) 

alaa  wincloaa(zfld_win); 

» 

if  (choica  **  12) 

( 

draw_fld  -  !draw_fld; 
if  (draw_fld) 

( 

minaiza (xmax/ 10, ymax/ 10) ; 
maxaiza (xmax-200, ymax-1 60) ; 
kaapaapact (xmax, ymax) ; 
fld_win  ■  winopan ("") ; 
wintit la ( fld_win_nama) : 
winconatrainta ( ) ; 
winaet (fld_win) ; 
zbuf far (TRUE); 
doublabuffer O  ; 

RGBmoda ( ) ; 
gconfigO; 
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shademodel (FLAT) ; 
miRodeCMVl  EWING); 


i 

els*  winclose ( f Id_win) ; 

I 

if  {{choice  >•  13)44(choice  <»  16)) 
transchoice  ■  choice  -  12; 
if  (tran3_choice  »“  4) 

I 

dxr  «  0; 
dxt  -  0; 
dxz  “  0; 
dyr  *  0; 
dyt  -0; 
dyz  •  0; 
tx  -  0; 
ty  •  0; 
t:  ■  0; 

trans_choice  «  1; 

1 

if  {choice  17) 


if 

(cur 

win 

seg  win) 

draw  seg  => 

!draw_seg 

if 

(cur 

win 

— 

cad_win) 

draw  rad  ■ 

idraw  rad 

if 

(cur 

win 

rat  win) 

draw  rat 

!draw_rat 

if 

(cur 

win 

con  win) 

draw  con  » 

idraw  con 

if 

(cur 

win 

— 

real_win)draw_real» 

draw  real 

if 

(cur 

win 

imag  win)draw  imag~ 

draw  imag 

if 

(cur 

win 

mag  win) 

draw  mag  - 

idtawmag 

if 

(cur 

win 

ph_win) 

draw  ph  »  idraw  ph; 

if 

(cur 

win 

»» 

eth_win) 

draw  eth  = 

!draw_eth 

if 

(cur 

win 

eph  win) 

draw  eph  » 

*draw_eph 

if 

(cur 

win 

-* 

zfld  win)draw  zfld- 

draw_2f Id 

if 

(cur 

win 

-* 

fid  win) 

draw  fid  « 

(draw  fid 

wine  lose (cur_win) ; 

) 

if  (choice  »»  18)  draw_axes  ■  'draw^axes; 
redraw  needed  •  TRUE; 

I 


) 


else 

( 


else  if  (dev  —  LEFTHOUSB) 

( 

if  (value) 

/•  read  initial  mouse  x,y,  which  will  be 

*  next  2  events  because  tied  to  leftmouse 

I 

qread((x) ; 
qreed(4y) ; 

/*  now  listen  to  all  mouse  x,y 
*  as  long  as  mouse  down  */ 
qdevice(H(XlSEX) ; 
qdevice(HOUSEY) ; 

I 

else 

/*  as  soon  as  mouse  released, 

•  stop  listening  to  mouse  x,y  •/ 

( 

unqdevice(HOUSEX) ; 
unqdevice (HOUSEY) ; 

I 

I 

if  (dev  *=  MIDDLEMOUSE  44  value  **  0) 

/•  do  a  pick  correlation  in  current  windcs-/  */ 
winset (cuc_win) ; 


data  choice  * 

»; 

if 

(cur 

win 

3eg_win) 

data^^choice  -  1; 

if 

(cur 

win  ■* 

rad  win) 

data  choice  *  2; 

it 

(cur 

win  *•= 

rat_win) 

data  choice  =  3; 

if 

(cur 

win  «« 

con  win) 

data  choice  =4; 

if 

(cur 

win  «* 

real  win! 

data  choice  ■ 

if 

(cur 

win 

imagwin) 

data  choice  »  6; 

if 

(cur 

win  == 

mag  win) 

data  choice  -  7; 

if 

(cur 

win 

ph  wir.)  data_choice  *  0; 

if 

(cur 

win 

eth  win) 

data^choice  =  9; 

if 

(cur 

win 

eph  win) 

data  choice  •  10; 

if 

(cur 

win  *  = 

zfld  win) 

data  choice  «  11 

if 

(cur 

win 

fid  win) 

data  choice  •  12; 

if 

(data  choice  11} 

tlireshtt; 

if  (rthresh  *1)  rthresh  *  0; 

) 

ei5»»  if  (data  choi'~e  -  =  !.■■) 

( 

1 1  h t esh* ♦ ; 

if  (fthresh  •*  *>)  tthresh  *  0; 

I 

else  if  ((data  choice  >«  l)44(data  choice  <»  6)) 


pi.:k(hiiffer,BUFSIZE); 
picking  •  TRUE; 
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deltax  -  S; 
d«ltfty  *  5; 

picksize(d«ltax.  daltay) ; 
iblink  »  0; 
if  (data  choice  >  0) 

( 

drawscene (dxr , dxt , dxz , dy r  «dyt « dyz « 

data  choice,  trana^choice,  iblink, 
buffer,  hits,  picking,  cur  win); 

) 

picking  *  FALSE; 
hits  •  endpick (buffer); 
if  (hits  >  0) 

{ 

iblink  •  buffertll; 
if  (draw_data  »•  FALSE) 

( 

prefposition(xaiax/40,xaax*8/20,yeiax/40, 
yMX*13/40); 
data_win  -  winopen(**“) ; 
wintitle(data_win_na«e) ; 
winconstraintsO ; 
winset (data_win) ; 
ainglebuf fer () ; 

RGBmode ( ) ; 
gconfigO ; 
shademodel (FLAT) ; 
anode  (MVlEVriNG); 
drawdata  -  TRUE; 

I 

winset (data^win) ; 

c3i(dbluecol); 

clear  () ; 

c3i  (vd^itecol) ; 

ortho2(0.0,  5,0,  0.0,  4.0); 

printhits ( iblink) ; 

I 

else 

I 

if  (draw  data) 

( 

dcaw^data  ■  FA1*SE; 
winclose(data  win); 

1 

f 

\ 

redraw  needed  ■  TRUE; 

) 

else  if  (dev  HOUSBX) 

( 

oldx  •  x; 

X  *  value; 

if  (tran3__choice  ••1)  dxr  •  dxr  ♦  (x-oldx); 

if  (trans^choice  »»  2)  dxt  -  dxt  ♦  (x-oldx); 

if  (trans_choice  •-  3)  dxz  -  dxz  ♦  (x-oldx); 

redraw_needed  •  TRUE; 

I 

else  if  (dev  -«  MOUSEY) 

( 

oldy  •  y; 
y  •  value; 

if  (trans_choice  ■■  1)  dyr  ■  dyr  ♦  (y-oldy); 

if  (tran3_choice  2)  dyt  ■  dyt  ♦  (y-oldy); 

if  (trans_choice  ••  3)  dyz  ■  dyz  ♦  (y-oldy); 

redraw_needed  «  TRUE; 

) 

)  /*  end  while  exitflag  FALSE 

and  (qtest  or  redraw^needed  FALSE)  */ 
)  /•  end  while  exitflag  FALSE  •/ 

exit (0) ; 

!  /♦  end  nalnO  */ 


/•  printhits 

*  Prints  wire  paraanter  data  in  data  window 

V 

void  printhits (iblink) 
long  iblink; 

( 

int  i , j ; 
char  str[60]; 
i  •  iblink-1; 
j  -  jwirelil-l; 
cmovfO. 2, 3. 6, 0. 0) ; 

sprintf(str,  "Model  contains  lid  wires,  with  lid  segments.",  nwires,  nsegs); 

charstr (str) ; 

ciaov(0.2,3.35,0.0); 

sprintf(str,  "You  have  selected  wire  lid,  segment  lid",  jwiredl,  iblink); 
charstr (str) ; 
cnov (0.2, 3. 1,0.0); 

sprintf(str,  "  (this  is  segiaent  lid  of  lid  on  TAG  number  tld).", 
j  seg ( l ) , isegs 1 j ) , it ags ( j I ) ; 
charstr (str) ; 
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cmov (0 . 2, 2. 85,0.0); 

sprintf(3tr,  "Wire  endpoints  ere:  (l.2f,  *.2f, 

bgnwire ( j ] [0 | , bgnwire ( j ] [ 1 ] , bgnwire [ j ) (2) ) ; 
charstr (str) ; 
cmov (0 . 2 , 2. 6, 0. 0) ; 

sprine£(str,  "  lt.2f.  I.2f,  l,2f>.", 

endwirel j  j 10), endwire ( j 1  1 1 } , endwire ( j 1 (2) ) ; 
c her str (str) ; 
cmov (0. 2,2.35,0.0) ; 

3printf(str,  "Segment  endpoints  ere;  {t.2f,  l.2f, 

bgnnode [ i ] ( 0 1 , bgnnode ( i ] 1 1 ] , bgnnode (i ) (2) ) ; 
cherstr (str) ; 
cmov 10.2,2.1,0.0); 

sprintf(str,  "  {l.2f,  l.2f, 

endnode ( i ] ( 0 ) , endnode ( 1 ) ( 1 ] , endnode [1)12)); 
cherstr  (st r ) ; 
cmov (0.2, 1.85,0.0) ; 

3printf(3tr,  "Segment  Length  -  t.2f  meters",  segsdl); 
charstr  (st r ) ; 
cmov 1 0. 2, 1 . 6, 0. 0) ; 

sprintf(str,  "Segment  Radius  *  t.3f  meters",  rads(i)); 

charstr  (str) ; 

cmov (0.2,1.35,0.0); 

sprintf(str,  "Segment  to  Radius  Ratio  -  srrat(i|}; 

charstr  (st r ) ; 
cmov (0.2, 1,1, 0.0); 
if  (dated)  (3)  «  0) 

sprint£(str,  "Both  ends  of  segment  are  connected."); 
if  (dated)  (3)  »=  2) 

sprintf(3tr,  "One  end  of  segment  is  connected."); 
if  (data(ij (3)  =»  4 ) 

sprintf(str,  "Neither  end  of  segment  is  connected."); 
charstr (str) ; 
cmov (0.2,0.85,0.0); 

sprintf(3tr,  "Frequency  (MH;):  l.2f",freq); 

■"harstr(str)  ; 
cmov (0.2, 0.6, 0.0); 

sprintffstr,  "Current  on  segment  freal,imag);  (l.4£,  l.4f)", 
rcr  I  i ) ,  icr  ( i ) ) ; 

charstr (str) ; 
cm<jv  (0.2,0.35,0.0); 

sprintfistr,  "Current  on  segment  (mag, phase):  t.4f)", 

mcr ( i] , xdata [ i I ) ; 

charstr (str ) ; 


/•  hsv2tgb  converts  from  hsv  to  rgb  colors.  It  assumes  that 

*  5  and  V  both  equal  1.  This  can  be  easily  changed  in  the 

*  future.  Result  is  put  into  global  variable  rgbcol(3). 

•/ 

void  hsv2rgb(f loat  hue) 

( 

float  h, 5, V, £ , p, q, t ; 
iiii  j; 

s  =  1.0; 

V  *  1.0; 
h  *  hue/ 60; 
j  =*  h; 
f  »  h-j; 
p  -  V  •  ( I  -  s ) ; 
q  «  V* ( 1  -  (s*  f ) )  ; 
t  -  vMl-(s*  (1-f))); 
switch  (j) 

( 

case  0: 

rgbcol ( 0 ) =v*  255 ; 
rgb..olIlJ=t*255; 
rgbcol (2 1 =p*  255; 
brealr; 
case  1 : 

rgbcol ( 0 ) *q*  255 ; 
rgbcol ( 1 ) *v*  255; 
rgbcol(21=p*255; 
break; 
case  2: 

rgbcol 1 0) -p*  255; 
rgbcol ) 1 ) *v*  255; 
rgbcol(2)*t*255; 
break; 
case  3; 

rgbcol 

rgbcol | 1 ) =q*255; 
rgbcol ( 2 ) -V*  255 ; 
break; 

rase  4 : 

rgbcol ( 0 1 • 255; 
rgbcol ( J |-p*255; 
rgbcoU21  =v*255; 
break ; 
case  5: 

rgbc"'!  I  n  ]  -v*  255; 
rgbcol ( 1 ] -p*  255 ; 
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rqbcol (2) ■q*255; 
br««k; 


/♦  drawscene 

*  drawscene  performs  all  drawinq.  Scene  is  cleared  and 

*  entire  scene  is  redrawn  each  time  called.  House  input 

*  is  used  to  change  polarview  viewing  transformation. 

* 

*  The  position  of  the  cube  is  fixed. 


void 

drawscene (short  dxr«  short  dxt,  short  dxz« 
short  dyr,  short  dyt,  short  dyz, 
long  data_choice,  long  trans_choice,  long  iblinJc, 
short  bufferO,  long  hits.  Boolean  picking,  long  do_win) 
i 

short  azimuth,  incidence,  fovy; 
float  cdr,  azi,  inc,  hue,  xstep; 
float  valmax,  valmin; 
int  i; 

unsigned  precision; 

static  long  vdiitecol  ( )  *  |  25S,  255,  255 
static  long  blackcol(]  *  {  0,  0,  0  ); 
static  long  redcold  •  (  255,  0,  0 
static  long  orangecol(]  »  (  255,  63,  0  ); 

static  Icmg  yellowcolt]  «  |  255,  255,  0 

static  long  greencol(]  «  |  0,  255,  0 
static  long  cyancolll  ■  (  0,  255,  255  1; 
static  long  bluecolf]  •  I  0,  0,  255  ); 
static  long  newbluecol()  >  I  64,  0,  64 
static  long  violetcolO  =  {  255,  0,  255 

static  long  greycol  ( )  >■  I  192,  192,  192  ); 


static 

float 

vbox 1(3) 

-  1  0.1,  0.2, 

0  ) 

static 

float 

vbox2t31 

-  {  1.3,  0.2, 

3  1 

static 

float 

vbox3(31 

•  1  1.3,  1.1, 

3  i 

static 

float 

vboxi (3) 

-  (  0.1,  1.1, 

3  ) 

static 

float 

vbox5(31 

-  {  1.3,  0.5, 

3  1 

static 

float 

vbox6(3) 

•  (  0.1,  0,5, 

3  ) 

static 

float 

vbox7(3) 

-  (  1.6,  0,2, 

3  ) 

static 

float 

vbcx8 ( 3 ] 

-  (  1.6,  0,7, 

3  I 

static 

float 

vbox9 ( 3 1 

-  {  0.1,  0.7,  0  1 

static 

float 

vl  ( 3 )  • 

0,2,  0.9,  0  ) 

static 

float 

v2(3l  - 

0,4,  0.9.  0  I 

static 

float 

v3(3l  - 

0,2,  0,8,  0  ) 

static 

float 

v4(31  - 

0.4,  0,8,  0  ) 

static 

float 

v5l3)  - 

0.2,  0,7.  0  1 

static 

float 

v6(3)  ■ 

0.4,  0.7,  0  ) 

static 

float 

v7t3l  - 

0.2,  0.6,  0  ) 

static 

float 

v8(3l  - 

0,4,  0,6,  0  » 

static 

float 

v9(3|  - 

0.2,  0,5.  0  ) 

static 

float 

vl0|3)  - 

(  0.4,  0.5.  0 

; 

static 

float 

vlH3)  - 

1  0,2,  0.4.  0 

; 

static 

float 

vl2(3J  * 

o 

o 

o 

; 

static 

float 

V1313)  - 

(  0.2,  0,3,  0 

; 

static 

float 

V1413)  - 

1  0,4,  0.3,  0 

; 

char  text (601; 


winset(do  win); 
c3i (blackcol ) ; 

if  (picking  FALSE)  clear (}; 

cdr  •  M_PI/10OO.O; 
fovy  »  450  ♦  dyz; 
if  (fovy  <  10) 

fovy  •  10; 

else  if  (fovy  >  ITOO) 
fovy  -  1790; 

azimuth  ■  dxr  *  5; 
incidence  •  dyr  •  5; 

if  (trans  choice  —  2) 

( 

azi  ■  cdr*azimuth; 
inc  ■  cdr*incidence; 

tx  ■  (dxt* fcos (azi)  -  dyt* fcos (inc)*fsin(azi) )/l00; 
ty  »  (dxt* f sin  (azi )  ♦  dyt* fcos ( inc) ♦ fcos (azi) )/ 100; 
tz  ■  (dyt* f sin (inc) )/ 100; 

I 

perspective ( fovy,  5. 0/4.0,  1.0,  9.0); 

if  (picking  ••  FALSE)  czclear (0x000000,  zval); 

pushmatrixO ;  /*  save  ModelView  matrix  by  pushing  and  duplicating 
polarview(5.0,  azimuth,  incidence,  0); 

/*  viewing  commands  premultiply  the  ModelView  matrix  ♦/ 
c3i (whitecol) ; 
pushmatrixO ; 

translate (tx, ty , tz) ; 
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wRiodel  (datd_choice,  iblink): 
popmdtcixO  ; 
c3i (greycol) ; 

if  ((picking  —  FALSE)  (dr«w_axes  —  TRUE))  axesO; 
popmatrixO;  /•  restore  ModelView  matrix  by  popping  ♦/ 

/♦  change  projection  to  ortho2  to  display  text  in  2D 

*/ 

ortho2(0.0,  5.0.  0.0,  4.0); 
linewidth  (2); 
if  (data_choice  ««  4) 

I 

/♦  make  wire  connections  key  */ 

/* 

bgnclosedline ( ) ; 
v3f (vboxl ) ; 
v3f (vbox7) ; 
v3f (vbox0) ; 
v3f (vbox9) ; 
endclosedline ( } ; 

•/ 

spr int f (text , 

•VIRE  CONNECTIONS’*); 
cinov(0.2,  0.6,  0.0); 
charstr (text ) ; 
c3i (newbluecol ) ; 
bgnlineO  ; 
v3f (v9) ; 
v3f (vlO) ; 
endline () ; 
c3i (whitecol) * 
sprintf (te 

•’Both  ends  C”  ected"); 
cmov (0  ,  . ■  ,  0.0); 

char'  r  '  .w) ; 

c3i  f'*:««»ncol) ; 
h-  liineO  ; 
v3f  (vll); 
v3f (vl2) ; 
endline ( ) ; 
c3i (whitecol) ; 
sprintf (text. 

“One  end  connected**); 
cmov (0.5,  0.4,  0,0); 
charstr (text) ; 
c3i (orangecol) ; 
bgnline( ) ; 
v3f (vl3) ; 
v3f (vl4); 
endline  0 ; 
c3i (whitecol) ; 
sprintf (text , 

"Neither  end  connected*’); 
cmov (0.5,  0.3,  0.0); 
charstr (text) ; 

) 

else 

! 

if  ((data  choice  <  8)  I  I  (-lata  choice  **  lUIUdata  choice  ==  12)) 

I  _ 

precision  *  2; 
if  (data_choice  **  1) 

I 

valmin  «  segmin; 
valmax  *  segmax; 

1 

if  (data_choice  «■  2) 

I 

precision  »  3; 
valmin  *  radmin; 
valmax  ■  radmax; 

) 

if  (data_choice  ==  3) 

( 

valmin  srmin; 
valmax  >  srmax; 
precision  -  0; 

) 

if  (data_choice  =*  5) 

( 

valmin  =  rcrmin; 
valmax  »  rcrmax; 

I 

if  (data_choice  *>=  6) 

( 

valmin  ■  icrmin; 
valmax  =  icrmax; 

I 

if  (data_choice  **  7) 

( 

valmin  «  mcrmin; 
valmax  *  mccmax; 

I 
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if  (data_choice  «*  11) 

{ 

valntin  -  0. ; 
valmax  •  zfldmax; 

if  (data_choice  --  12J 

1 

valmin  *  0.: 
valmax  ■  fldnax; 

) 

if  (valmin  ••  valmax) 

( 

bgnclosedlineO ; 
v3f (vboxl) ; 
v3f (vbox2) ; 
v3f (vboxS) ; 
v3f (vbox6) ; 
endclosedline ( ] ; 

if  (data_choice  —  1)  aprintf (text, “SEGMENT  LENGTH,  METERS“); 
if  (data  choice  ■-  2)  sprintf (text, "WIRE  RADIUS,  METERS"); 
if  (data^choice  —  3)  aprintf (text, "SEGMENT  TO  RADIUS  RATIO"); 
if  (data  choice  —  5)  sprintf (text, "REAL  COMPONENT  OF  CURRENT"); 
if  (data_choice  —  6)  sprint f (text, "IMAG  COMPONENT  OF  CURRENT"); 
if  (data  choice  —  7)  sprlntf (text, "MAGNITUDE  OF  CURRENT"); 
if  (data  choice  11)  sprintf (text, "2-COMP  OF  NEAR  FIELD,  VOLTS/M") 
if  (data__choice  ■■  12)  sprintf  (text,  "TOTAL  NEAR  FIEID,  WLTS/M"); 
cmov (0.2,  0.4,  0.0); 
charstr (text) ;  ^ 

c3i (newbluecol) ; 
bgnlineO ; 
v3f (vl3) ; 
v3f (vl4); 
endline  () ; 
c3i (whitecol) ; 
sprintf (text, 

"16. *f",  precision,  valmin); 
cmov (O.S,  0.3,  0.0); 
charstr(text) ; 

) 

else 

( 

bgnclosedline ( ) ; 
v3f (vboxl) ; 
v3f (vbox2) ; 
v3f (vbox3) ; 
v3f (vboxl ) ; 
endclosedline ( ) ; 

if  (data  choice  1)  sprintf  (text,  "SE01BNT  LENGTH,  METERS"); 

if  (data^choice  2)  sprint  £ (text , "WIRE  RADIUS,  METERS"); 

if  (data'^choice  *»  3)  sprintf  (text.  "SEGMENT  TO  RADIUS  RATIO"); 

if  (data  choice  —  5)  sprintf (text. "REAL  COMPONENT  OF  CURRENT"); 

if  (data'choice  ««  6)  sprintf (text. "IMAG  COMPONENT  OF  CURRENT"); 

if  (data  choice  ■*  "))  sprintf  (text.  "MAGNITUDE  OF  CURRR^T"); 

if  <data_choice  »»  11)  sprintf (text, "2-CCWP  OF  NEAR  FIELD,  VOLTS/M"); 

if  (data'choice  ■»  12)  sprintf (text, "TOTAL  NEAR  FIELD,  VOLTS/M"); 

cnov{0.2,  1.0,  0.0); 

charstr (text) ; 

sprintf (text, 

"16. *f  to  16. ♦f",  precision,  valmin,  precision, 
valmin* (val«ax-valmin)/7) ; 

cmov (0.5,  0.9,  0.0); 
charstr (text ) ; 
sprintf (text, 

"16. *f  to  16. 'f",  precision,  valmin* (valmax-valmin) /T, 
precision,  valmin+2' (valmax-valmin) /7) ; 
cmov(0.5,  0.8,  0.0); 
charstr(text)/ 
sprintf (text, 

"16. 'f  to  16.*f",  precision,  valmin*2*  (valiBax-val«in)/7, 
precision,  valmin* 3* (valmax-valmin) /7) ; 
cmov(0.5,  0.7,  0.0); 
charstr (text) ; 
sprintf (text, 

"16. *f  to  16. ♦f",  precision,  val»in*3* (valmax-valmin) /7, 
precision,  valmin*!' (va Imax- valmin ) /7) ; 
cmov (O.S,  0.6,  0.0); 
charstr (text) ; 
sprintf (text , 

"16. *f  to  16. 'f",  precision,  valmin*!* (valmax~valrain)/7, 
precision,  val»in*S* (valmax-valmin)/7>; 

cmov (0.5,  0.5,  0.0); 

charstr (text) ; 
sprintf (text, 

"16. ♦f  to  16.*£",  precision,  valmin*5* (valmax-valmin)/7, 
precision,  valmin+6* (valmax-valmin)/7) ; 
cmov (0 . 5,  0. ! ,  0.0); 
charstr (text ) ; 
sprintf (text, 

"16. *f  to  16.*f",  precision,  vaimin*6* (valmax-valmin)/7, 
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precision,  vaijnin+  {vaJmax-valmi  n) )  ; 
cmQv(0.S,  0.3,  0.0); 
charstr (text) ; 
if  (data_choice  <  8) 

( 

c3i (newbluecol) ; 
bgnline  ( ) ; 
v3f (vl); 
v3f  (v2)  ; 
endline ( } ; 
c3i (cyancol) ; 
bgnline  ( ) ; 
v3f (v3); 
v3f(v4); 
endline  () ; 
c3i (greencol ) ; 
bgnline ( ) ; 
v3f (v5) ; 
v3f (v6) ; 
endline ( ) ; 
c3i (yellowcol ) ; 
bgnline ( ) ; 
v3f(v7); 
v3f (v8) ; 
endline  ( ) ; 
c3i (orangecol) ; 
bgnline ( ) ; 
v3f {v9) ; 
v3f (vlO) ; 
endline  ( ) ; 
c3i ( redcol) ; 
bgnlineO  ; 
v3f (vll) ; 
v3f (vl2) ; 
endline  ( ) ; 
c3i  (v^itecol ) ; 
bgnlineO  ; 
v3f {vl3> ; 
v3f (vl4 ) ; 
endline ( ) ; 
c3i (whitecol) ; 

) 

else 

/♦  near  fields  color  key  •/ 

I 


c3i (newbluecol) 
rectf(0.2,  0.9, 

0.4, 

0.96) 

c3i (cyancol) ; 
rectf(0.2,  0.8, 

0.4, 

0.86) 

c3i (greencol) ; 
rectf(0.2,  0.7, 

0.4, 

0.76) 

c3i (yellowcol) ; 
rectf(0.2,  0.6, 

0.4, 

0.66) 

c3i (orangecol ) ; 
rectf(0.2,  0.6, 

0.4, 

0.66) 

<.3i  (redcol) ; 
rectf(0.2,  0.4, 

0.4, 

0.46) 

c3i (whitecol) ; 
rect f (0.2,  0.3, 

0.4, 

0.36) 

c3i (whitecol ) ; 

) 

1 

} 

if  (data_choice  *=  11) 

( 

if  (zthresh  <=  6) 

I 

sprintf  (text, ’’THRESHOLD  CUT  AT  16. 2f  VOLTS/M",  zthresh*valmax/7. ) 
1 

else  if  (zthresh  ==  7) 

( 

sprintf (text,  "THRESHOLD  CUT  AT  HERP  (16. 2f  VOLTS/M)",  herp); 

} 

else  if  (zthresh  ==  8) 

I 

sprintf  (text, "THRESHOLD  CUT  AT  HERO  (46. 2f  VOLTS/M)’’,  hero)  ; 

) 

cmov(0.2,  3.86,  0.0); 
charstr (text) ; 

sprintf (text, "CLICK  ON  MIDDLE  MOUSE  TO  CHANGE  THRESHOLD”) ; 
cinov(0.2,  3.7,  0.0); 
charstr (text) ; 

spr int f (text , "FREQUENCY  =  4. If  MHz",  freq); 
ctnov (0. 2,  3-65,  0.0); 
charstr (text) ; 

) 

if  (data_choice  12) 

\ 

if  (tthresh  <*  6) 

( 

sprintf (text, "THRESHOLD  aJT  AT  46. 2f  VOLTS/M",  tthre3h*valnax/7. ) 
) 

else  if  (tthresh  ==  7) 
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I 

sprintf (text, "THRESHOLD  CUT  AT  HERP  /%6.2f  VOLTS/M)",  herp) ; 

I 

else  if  (tthresh  »•  8) 

( 

sprintf (text, "THRESHOLD  CUT  AT  HERO  (16. 2f  VOLTS/M>",  hero); 

} 

cmov(0.2,  3.85,  0.0); 
charstr (text) ; 

sprintf (text, "CLICK  ON  MIDDLE  MOUSE  TO  CHANGE  THRESHOLD"); 
cmov(0.2,  3.7,  0.0); 
charstr(text) ; 

sprintf (text, "FREQUENCY  -  l.lf  MHz",  f req) ; 
cinov(0.2,  3,55,  0.0); 
charstr (text) ; 

) 

if  (data_^choice  <®  8) 

( 

sprintf (text, "USE  MIDDLE  MOUSE  TO  SELECT  A  WIRE  SEGMENT"); 
cmov{0.2,  3.85,  0.0); 
charstr (text) ; 

} 

if  ( (data_choice  >*  5) 4« (data_choice  <»  8)) 

( 

sprintf (text, "FREQUENCY  =  t.lf  MHz",  freq); 
cmov (0.2,  3.7,  0.0); 
charstr (text) ; 

) 

if  ( (data_choice  *=  8) I  I (data_choice  =*  9) I  1 (data_choice  «»  10)) 
( 

if  (data  chov'-e  =*  9) 

( 

sprintf(text, "MAX  MAGNITUDE:  t6.2f",  ethmax) ; 
cmov(0.2,  3.85,  0.0); 
charstr (text) ; 

sprintf (text, "FREQUENCY  *  l.lf  MHz",  freq); 

C7nov(0.2,  3.7,  0,0); 
charstr (text) ; 

I 

if  (data_choice  *»  10) 

I 

sprintf (text, "MAX  MAGNITUDE:  16, 2f",  ephmax) ; 
cmov(0,2,  3,85,  0.0); 
charstr (text) ; 

sprintf (text, "FREQUENCY  »  l.lf  MHz",  freq); 
cmov(0,2,  3.7,  O.OJ; 
charstr (text ) ; 

I 

sprintf (text, "PHASE") ; 
cmov(0.2,  1.0,  0,0); 
charstr (text) ; 
for  (i*0;  i<30;  i*+) 
i 

hue  •  360.  -  i*6.-3.; 

hsv2r9b(hue) ; 

c3i (rgbcoi) ; 

xstep  »  0.2  +  i  *  .03; 

rectf(xst€p,  0.7,  xstep>.03,  0.85); 

) 

for  {i»0;  i<30;  !♦+) 

{ 

hue  =  180 .  -  i* 6, - 3 . ; 

hsv2rqb(hue) ; 

c3i (rgbcoi) ; 

xstep  »  0.2  +  i  •  .03; 

rect f (xstep, 0.5, xstep* . 03, 0. 65) ; 

I 

c3i (whitecoi) ; 
sprintf (text, "-180") ; 
croov(0. 2, 0. 9, 0.0); 
charstr (text) ; 
sprintf (text, "0") ; 
cmov (0.2, 0.4, 0.0); 
charstr (text) ; 
sprintf (text, "0") ; 
cmov (1.06,0.9,0.0) ; 
charstr (text) ; 
sprintf (text, "180") ; 
cmov (0 . 98,0.4,0.0); 
charstr (text ) ; 

) 

I 

c3i (whitecoi ) ; 
if  (trans_choice  -»  1) 

I 

sprintf (text, 

"LEFTMOUSE  ROTATES:  (I5d, I5d, I4d, " 

"I5.2f,l5.2f,t5.2f)", 

azimuth,  incidence,  fovy,  tx,  ty,  tz); 

) 

if  (ttans  choice  *■  2) 

( 

sprintf (text, 
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(FILE*)0) 


if((etpfil«  ■  fop«n(etp  « 

( 

printf ("\n\t\t  File  Is  could  not  be  opened.  Nn**,  etp_file); 
exit  ( ) ; 

) 

i  •  0; 
iphis  -  0; 
ithetas  -  0; 

fscanf (et»f ile,  "Is*’,  line)' 
fscanf (etnf lie,  line): 

fscanf (etmfiie,  “If  ,  If  ,  If  ,  If  ,  If  ,  If”, 
cminnag,  cethmax,  4a,  4b,  4c,  4d); 
fscanf (etnfile,  ”lf  ,  If  ,  If”,  4ffiag.  4phi,  4theta); 
oldtheta  •  theta; 
oldphi  ■  phi; 

while  ((mag  >  -1.23)||(mag  <  -1.24)l|(phi  >  -1.23)ll(phi  <  -1.24) 

I  I (theta  >  -1.23) I  I (theta  <  -1.24>) 

i 

if  (theta  >  oldtheta)  ithetas44; 
if  (phi  >  oldphi)  iphis++; 
if  (theta  <  oldtheta)  ithetaa  »  0; 
if  (phi  <  oldphi)  iphis  ®  0; 

ethdata  (iphis)  (i thetas)  (0)  *  2*iiiag*  fcos  (cdr* theta)  *  fsin(cdr*phi}/ethAax; 
ethdata (iphis) (ithetasj (1)  *  2*mag* fsin(cdr*theta) • f3in(cdr*phi) /ethmax; 
ethdata ( iphis] I ithetas] [2]  -  2*mag* fcos (cdr*phi)/ethmax; 
i++; 

oldtheta  ■  theta; 
oldphi  «  phi; 

fscanf (etfflfile,  "If  ,  If  ,  If",  4niag,  4phi,  4theta); 

) 

fclose (etmfiie) ; 

npts  »  i; 

nphis  *  iphis  *  1; 

nthetas  «  ithetas  ♦  1; 

fscanf (etpfile,  "Is",  line); 

fscanf (etpfile,  "Is",  line); 

fscanf (etpfile,  "If  ,  If  ,  If  ,  If  ,  If  ,  If",  4a,  4b,  4C,  4d,  4e,  4f); 
for  (i»0;  i  <  nthetas;  !+♦) 

{ 

for  (3»0;  j  <  nphis;.  j  +  +  ) 

/ 

fscanf (etpfile,  "If  ,  If  ,  If",  4phase,  4phi,  atheta); 

/*  put  phase  between  0  and  360  */ 

while  (phase  <  0.0)  phase  »  phase  *  360; 
while  (phase  >  360.0)  phase  »  phase  -  360; 
ethphase( j] (i)  «  phase; 

) 

) 

fclose (etpfile) ; 

eth  initialized  «  TRUE; 


/*  initialize  the  E-phi  component  data  ♦/ 
static  void  initeph (void) 

( 

float  a,  b,  c,  d,  e,  f; 

float  mag,  phase,  phi,  theta,  minmag,  oldphi,  oldtheta; 
float  cdr; 

int  i,  j,  iphis,  ithetas; 
char  line['70]; 
cdr  •  M_PI/180.0; 

if((epmfile  *  fopen (epm_file, "r") >  =«  (FILE*)0) 

{ 

printf ("\n\t\t  File  Is  could  not  be  opened.  \n",  epin_file); 
exit  O  ; 

) 

if((eppfile  »  fopen (epp_file, "r") )  (FILE*)0) 

{ 

printf ("\n\t\t  File  Is  could  not  be  opened.  \n",  epp_file); 
exit  0 ; 

) 

i  -  0; 
iphis  -  0; 
ithetas  ■  0; 

fscanf (epmfile,  "Is",  line); 
fscanf (epeifile,  "Is",  line); 

fscanf  (epmfile,  "If  ,  If  ,  If  ,  If  ,  If  ,  If", 

4minmag,  4ephmax,  4a,  4b,  4c,  4d): 
fscanf (epmfile,  "If  ,  If  ,  If",  imag,  4phi,  4theta); 
oldtheta  -  theta; 
oldphi  “  phi; 

while  ((mag  >  -1.23)ll(mag  <  -1.24)n(phi  >  -1.23))Mphi  <  -1.24) 

IMtheta  >  - 1 . 23)  M  (theta  <  -1.24)) 

( 

if  (theta  >  oldtheta)  ithetas**; 
if  (phi  >  oldphi)  iphis**; 
if  (theta  <  oldtheta)  ithetas  »  0; 
if  (phi  <  oldphi)  iphis  •  0; 

ephdata (iphis) (ithetas) (0)  “  2*mag* fcos (cdr*theta) * fsin(cdr*phi) /ephmax; 
ephdata (iphis ) ( Ithetas 1 ( 1 ]  *  2*mag* fsin (cdr* theta) * f sin (cdr^phi) /ephmax; 
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»phdata(iphis) lithetasl 121  *  2*mag* fcos (cdr*phi) /ephnax; 

i++; 

oldtheta  •  theta; 
oldphi  •  phi; 

fscanf (epmfile,  "%f  ,  ,  »£",  tmag,  4phi,  4theta); 

I 

f close  (epiBfile) ; 
npts  •  i; 

nphis  »  iphis  ♦  1; 
nthetas  •  ithetas  +  1; 
fscanf (eppfile,  "Is",  line); 
fscanf (eppf iie,  "Is",  line); 

fscanf(eppfile.  "«f  ,  If  ,  If  ,  »f  .  «f  .  If".  “f' 

for  {i*0;  i  <  nthetas;  i*+) 

I 

for  (j-0;  j  <  nphis;  j++) 

fscanf (eppfiie,  "If  ,  If  ,  If**,  4phas€,  4phi,  4theta); 

/*  put  phase  between  0  and  360  */ 

while  (phase  <  0.0)  phase  =  phase  *  360; 
while  (phase  >  360.0)  phase  =  phase  -  360; 
ephphasel jl (il  -  phase; 

1 

) 

fclose(eppfile) ; 
eph_initialized  »  TRUE; 


4e,  4f); 


/♦  initialize  the  near  field  data 
static  void  initf Id (void) 

1 

float  X,  y,  z,  zfld,  fid; 
int  i,  j; 

if((fldfile  =»  fQpen(fld^file,"r’‘))  »»  (FILE*)0) 

*  printf ("\n\t\t  File  Is  could  not  be  opened.  \n",  fld^file); 
exitO  ; 

1 

zfldnax  -  0.; 
fldnax  *  0.; 

fscanf (fldfile,  "Id",  4nflds); 

for  {i»0;  i  <  nflds;  i+*) 

*  fscanf (fldfile,  "If  If  If  If  If",  4x,  4y,  *z,  ftzfld,  4fld); 
if  (zfld  >  zfldmax)  zfldmax  »  zfld; 

if  (fid  >  fldmax)  fldmax  •  fid; 
fldpntsdl  (01  *  tx-xshift) ‘mscale; 
fldpntsli) (1)  »  (y-yshift)*mscale; 
fldpntsUl  (21  -  {2-zshift)*mscaU; 
zflddatali)  -  zfldMOO.; 
flddatadl  -  fldMOO./ 


I 


1 

for  (i«0;  i  <  nflds;  !+♦) 
( 


tfldpntsdl  (0) 
tfldpntsdl  (I  j 
tfldpntsli) (2] 
zfldpntsd]  (OJ 
zf Idpnts ( i ] ( 1 ) 
zf Idpnts (i ] [2 1 
zspheres(i] (0) 
zspheres(i] (1] 

zspheres (il (21 
tspheresfi] (0] 
tspheresUJ  P) 
tspheres(i) (21 
zspheres ( i] [3] 
tspheres ( i ) ( 3 ) 


fldpnts ( i! (01 ; 
f Idpnts f i i  ( If ; 

fldpnts (i| (2 l+flddata(i)/ fldmax/ 100. ♦mscale; 
fldpntsli) |0); 
fldpnts ( il 1 1 1 ; 

fldpntsli  I  (21  ♦zflddatali) /zfldmax/ 100. *iB5cale; 

fldpntslil (01; 

fldpntsUl  (Il ; 

fldpnts ( i ) (2) ; 

f Idpnts (i 1 (01 ; 

fldpntsli) il); 

fldpntsUl  (2) ; 

Zflddata { i 1 /zf ldmax/4 00. ‘mscale; 
f Iddata (i 1 / f ldmax/400. ‘mscale; 


) 

fclose(fldfiie) ; 
£ld_initialized  »  TRUE; 


/*  initialize  the  models  vertices  •/ 
static  void  initmodel (void) 

I 

float  a,  b,  c,  d,  e,  f,  xpcr; 

int  i,  itag,  iseg,  Irad,  isr,  icon,  ircr,  iicr,  imct,  iw,  is; 


fscanf (geofi le,  "If",  4freq); 

if  (freq  <  3.)  herp  =  632. S; 

else  if  (freq  <  30.)  herp  »  1097. /freq; 

else  if  (freq  <  300.)  herp  ^  63.26; 

else  if  (freq  <  IbOO.}  herp  =  3.65  •  sqrtifjeq); 

else  herp  ®  141,4; 

if  (freq  <  1.)  hero  =  100. /freq; 

else  if  (freq  <  3.7)  hero  *  100 . /pow( f req, 3. ) ; 

else  if  (freq  <  10.)  hero  »  2.; 

else  hero  •  0.5*pow{freq,0.6); 

fscanf (qeofile,  "Id  Id",  4nwires,  4n5eg3); 
fscanf (geofile,  "If  If",  4mxmin,  4mxmaH); 
fscanf (geofile,  "If  If",  4mymin,  4mymax) ; 
fscanf (geofile,  "If  If",  4mzmin,  4mzmax); 
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fscanf (9eofile,  “tf  tf”,  exshift,  cyshift); 

fscanf (geof ila,  "tf  tf**,  czshift,  emscale); 

fscanf (geofile,  '*lf  tf",  csegmin,  esegmax); 

fscanf (geofile,  "tf  If**,  era^in,  eradmax); 

fscanf (geofile,  "%f  If",  esrmin,  esrmax); 
fscanf (geofile,  "If  If",  erermin,  erernax); 

fscanf (geofile,  "If  If",  ciernin,  ciernax); 

fscanf  (geofile,  "If  If",  (rnermin,  cinermax); 

for  (i  ■  0;  i  <  nwires;  i++) 

( 

fscanf (geofile,  "Id  Id  If  If  If  If  If  If",  citag,  eiseg,  ea, 
<b,  cc,  cd,  (e,  «f); 
itagsli]  »  itag: 
isegs[i]  «  iseg; 
bgnwire(i) (0]  »  a; 
bgnwire(i) [II  »  b; 
bgnwi  re ( i ] ( 2 ]  -  c ; 
endwire(il (0)  «  d; 
endwire{i] ( 1 ]  »  e; 
endwire(i) [2)  -  f; 

] 

for  (i  •  0;  i  <  nsegs;  i++) 

( 

fscanf (geofile,  "Id  Id  Id  Id  Id  Id  Id  Id  Id  If",  eis,  tiw, 
ciseg,  eirad,  cisr,  eicon, 

Circr,  ciicr,  eimer,  exper); 
jseg{i)  -  is; 
jwirelil  -  iw; 
data [i] (0)  «  iseg; 
data [i] ( 1 )  »  irad; 
data(l] [2]  •  isr; 
data(i] (31  -  icon; 
datalil (41  «  ircr; 
data(i]  [5)  =»  iicr; 
data (i) (6)  »  imer; 
xdata[i]  *  xper; 

fscanf (geofile,  "If  If  If  If  If  If",  ca,  eb,  (c,  ed,  ee,  ef); 

segs(i]  "  a; 

rad3(i)  -  b; 

srratti]  «  c; 

rcr (i]  •  d; 

icr(i]  >  e; 

nc r  ( i )  >  it 

fscanf  (geofile,  "If  If  If  If  If  If",  ea,  tb,  *c,  4d,  4e,  *f): 
bgnnodedi  (0]  •  a; 
bgnnode(i) (II  *  b; 
bgnnode ( i ) ( 2 ]  •  c ; 
endnode(i) (0|  ■  d; 
endnode[i] (1)  »  e; 
endnode ( i J ( 2 J  ■  f; 
bgndatadl  (0)  ■  (a-X3hift)*m3cale; 
bgndata(i] (I]  *  (b~yshift) *mscale; 
bgndata(il (2)  “  (c-zshift) ‘mscale; 
enddata [ i ] [0]  •  (d>xshift) *mscale; 
enddata (i ] ( 1 )  *  (e-yshift) ♦mscale; 
enddata (ij (2)  *  (f-zshift) ♦mscale; 

I 

fclose (geofile) ; 

model  initialized  >  TRUE; 


/♦  wireframe  NEC  model  */ 

void  wmodeKlong  data_choice,  long  iblink) 

( 

int  i,  j; 

float  h,s,v, f,p,q,t,val; 
static  float  originlJ  »  i  0,  0,  0  1; 
static  long  redcoK)  ■  (  255,  0,  0  J; 
static  long  orangecol()  *  (  255,  63,  0  ); 
static  long  yellowcoll)  •  (  255,  255,  0  1; 
static  long  greencold  ■  (  0,  255,  0  ); 
static  long  cyancold  •  |  0,  255,  255  1; 

static  long  bluecold  *  (  0,  0,  255  ); 

static  long  newbluecold  "  (  64,  0,  64  |; 
static  long  violetcolO  -  (  255,  0,  255  1; 

static  long  greycol ()  -  {  192,  192,  192  1; 

static  long  purplcold  -  (  255,  0,  176  J; 
static  long  whitecold  »  (  255,  255,  255  ); 
linewidth(2) ; 

if  (!  nodel^initialized) 
initmodel ( ) ; 

if  ((!  eth_initialized) 44 (data_choice  9)) 

initethO ; 

if  ((!  eph_initialized)44  (data^choice  10)) 

initeph( ) ; 

if  ((!  fld_initialized) 44 ( (data_choice  —  1 1 ) I  I (data  choice  12))) 
initfldO  ; 

/♦  draw  segments  ♦/ 
if  (data_choice  <  9) 
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for  U  ■  0;  i  <  ns«gs;  i++) 

{ 

if  (data_choice  <  0) 

if  (data(i)  (data_^choic«--l}  «®  0> 
c3i (newbluecol) ; 

if  (data(ij [data_choice-J )  =*  1) 
c3i (cyancol) ; 

alse  if  (data ( i} (d«ta_choic«-l 1  **  2) 
c3i (greencol) ; 

else  if  (data (ij (data_choice-l 1  **  3) 
c3i (yellowcol) ; 

else  if  (data ( i) (data_choice- 1 )  =-  4) 
c3i (orangecol) ; 

else  if  (datali) (data_choice-ll  »=  b) 
c3i (redcol ) ; 

else  if  (data t i } (datd_choice- 1 1  ==  6) 
c3i  (vdiitecol) ; 

J 

if  (data_choice  8) 

( 

hsv2rgb(xdata { i ] ) ; 
c3i (rgbcolj  ; 

1 

loadnamed*  1 ) ; 
if  (i  !■  ibiink-1) 

( 

bqnlineO; 

v3f (bgndata ( i ] ) ; 
v3f (enddata [1 } ) ; 
endline () ; 

I 

) 

» 

if  (data_choice  «*  9) 

( 

for  (i  «  0;  i  <  nphis-l;  i*+) 

\ 

bgnqstcipO  ; 

for  (j  ■  0;  j  <  nthetas;  j**) 

( 

hsv2rgb(ethphase( i I ( j I ) ; 

c3i (rqbcol) ; 

v3f (ethdata ( i I ( j  1  ) ; 

hsv2rgb(ethphase( i^l 1 ( j  I ) ; 

c3i (rgbcol) : 

v3f  (ethdata  ( i*l  j  ( j  1 )  -• 

( 

endqstripO  ; 

) 

I 

if  (data  choice  10) 

1 

for  (i  -  0;  i  <  nphis  1;  i**) 

I 

bgnqst rip ( ) ; 

for  (j  «  0;  i  <  nthetas;  3*>) 

1 

hsv2rgb (ephphase ( i I ( i  1 ) ; 

c3i (rgbcol) ; 

v3f (ephdata ( i | ( j I ) ; 

hsv2rgb (ephphase ( i *  1 ] [ j ) ) ; 

c3i (rgbcol; ; 

v3f (ephdata | i*l ) ( 3  1  )  ; 

I 

endqstripO  ; 

) 

I 

if  (data_choice  11) 

I 

c3i  (greycol) ; 

for  (i  -  0;  i  <  nsegs;  !*♦) 

{ 

bgnline  ( ) ; 

v3f (bgndata (i ) ) ; 
v3f (enddata  [  i ) ) ; 
endline  ( ) ; 

I 

for  (i  “0;  i  <  nflds;  !♦♦) 

( 

val  ■  z f Iddata ( i 1 / z f ldmax*7. / 100. ; 
if  (((val  >«  zthresh) *4 (zthresh  <  7))j| 

( (zf  Iddata  (  i  )  >-  herp*  1 00 . )  44  (zthresh  -• 

( (zf Iddata ( i I  >■  hero* 100. )44 (zthresh  *=  8))) 
I 

if  (val  <  1 . ) 

I 

c3i (newbluecol ) ; 
bgnpoint ( ) ; 

v3f (f Idpnts ( i 1 ) ; 
endpoint ( ) ; 

) 
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else 

{ 

initx  »  f Idpnts (i ) (0)  -  mscale/2.; 
inity  =  f  Idpnts  (  i  )  1 1 )  -  Riscale/2.; 
initz  •  fldpnts (i) (2)  -  m3cale/2.; 
if  (val  <  2.) 

c3i (cyancol) ; 
num  •  2: 

» 

else  if  (val  <  3.) 

I 

c3i (greencol) ; 
num  -3; 

) 

else  if  (val  <  4.) 

I 

c3i  (yellowcol) ; 
num  •  5; 

else  if  (val  <  5.) 

I 

c3i (orangecol > ;  • 

num  »  7; 

else  if  (val  <  6. ) 

I 

c3i (redcol) ; 
num  *  9 ; 

I 

else 

{ 

c3i  (vdiitecol) ; 
num  *  11; 

I 

step  *  (mscale/2. )/( (num-1 .) /2. ) ; 
for  (inx  *  0;  inx  <  num;  inx*+) 

( 

fpntIO)  «  initx  ♦  inx  *  step; 
for  (jnx  *  0;  jnx  <  num;  jnx++) 

( 

fpnt(ll  ■  inity  ♦  jnx  *  step; 
bgnpoint  O ; 

for  (knx  »  0;  knx  <  num;  knx**) 

( 

fpnt(2)  ■  initz  *  knx  ♦  step; 
v3f (fpnt) ; 

I 

endpoint (} ; 

) 

1 

) 

) 

} 

linewidth (2} ; 

I 

if  (data  choice  ««  12) 

( 

c3i (greycol) / 

for  (i  •  0;  i  <  nsegs;  i+*) 

( 

bgnline (} ; 

v3f (bgndata ( i} ) ; 
v3f (enddata ( 1 ))  ; 
endline ( ) ; 

) 

for  (i  ■  0;  i  <  nflds;  i++) 

( 

val  ■  f  iddata  f  i  J/ f  lclmax*7. /lOO.  ; 
if  (val  >»  6.)  c3i (whitecol) ; 
else  if  (val  >»  5.)  c3i (redcol); 
else  if  (val  >•  4 . )  c3i (orangecol) ; 
else  if  (val  >»  3.)  c3i (yellowcol) ; 
else  if  (val  >■  2.)  c3i (greencol) ; 

else  if  (val  >»  1.)  c3i (cyancol) ;  * 

else  c3i (newbluecol) ; 

if  (((val  >•  tthresh) *« (tthresh  <  7))|| 

((flddatalil  >-  herp* 100 . ) «« (tthresh  7))(| 

((flddata(i)  >»  hero* 100 . ) a* (tthresh  -*  8))) 


bgnline ( ) ; 

v3f (fldpnts (i) ) ; 
v3f (t fldpnts i i I ) ; 
endline ( ) ; 

if  (val  <  1.) 

( 

bgnpoint ( ) ; 

v3f (fldpntsli) ) ; 
endpoint ( ) ; 

) 

else 

sphdraw(tspheres li} ) ; 
I 

) 

linevidth(2) ; 

) 


I 
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Appendix  C 

SOURCE  CODE  FOR  view_mom 

view_mo«i.c 

view  mo*#  displays  a  ol  moments  matrix. 


•  ini'  1  ud**  /  g  1 .  h> 
•inzlude  <gl/device . h > 
i  in-:  iude  <'stdio.  h  ' 
•include  <mat h. h> 
•define  BUFSIZE  SO 


/•  function  prototypes 

*/ 

void  init  la  I  ize  ( V  ..id  >  ; 
void  main(void); 

void  pc  int  hit  3  t  long  iblink,  sh.'r*;  f.«,  sh-rt  •olumn); 

Void  hsv2 rgb { f loat  h); 

void  drawscene ( short  dxr,  short  dxt .  short  dxz,  short  dxs, 
short  dyr,  short  dyt ,  short  dyz,  short  dys, 

long  datachoice,  long  t ranschoioe,  long  display  node,  long  iblink, 
long  do  win) ; 
void  axes (void); 
void  ini treat c IX (void) ; 

void  momreat  r  i  X  ( 1  onq  data  choice,  Ir-ng  display  mode,  long  iblink); 

FILE  i  1  e; 

char  f i lename 125]; 

chat  iBom_f  i  le  ( 2b  1 ; 

static  char  mom[j  *  '’.mom”; 

int  nnodesx, nnodesy; 

float  mommag [ 1000 j  [ 1 000 1 ,  momphase [ 1 OOO ) 1 1000 1  ; 
float  datallOOO) llOOn]  (I); 
float  zscale,  mommax.  miTfiimin; 
long  tgbcolOl; 

Boolean  picking  rcM  -x  FALSE; 

Boolean  picking_col  «  FALSE; 

static  float  tx,  ty,  tz; 

/*  window  id’s  •/ 

long  mom_win,  ph^win,  cuc_win,  data_win; 
static  char  reo(B_win_name i !  »  ”MoM  Matrix"; 
static  char  ph_win_name ( )  •  "MoM  Matrix  with  Phase"; 
static  char  data_win_name ( )  *  "Pick  Correlation  Info”; 
long  xreax,  ymax,  zval; 


/ •  init  ia I  ize 

♦  The  initialize  subroutine  positions  the  window  and  specifies 

•  its  future  const  t  aint  .  Graphics  ‘-onf  igurat  ions  are  set  and 

*  the  event  queue  is  ini t ia 1  iced. 

•/ 

void 

init ill ize (void) 

( 

float  mreal,  reireaq,  ma-j,  phi.ce,  --.Ir; 
int  i.  j; 

printf("  Enter  file  nam*-:  "); 

3Canf("ls",  filenime); 

St  rcpy  (mom  f  i  le,  f  i  1  e_nani*'j  ; 

St r cat  (mocn_file,mom)  ; 

if|{momfile  »  fopen(mom  file,'*r"))  (FILE*)0) 

I 

printf  ('’\n\t\t  File  Is  c<*iil-l  not  be  opened.  \n",  mom_file); 
exit  ()  ; 

) 

cdr  '  M_PI/180.0; 
mtunmax  ^  0; 
moiwnin  =  lOOOOO; 

fscanf  (momf  i  le,  "Id  Id",  <nn*»de5x,  4nnod*?5y)  ; 
for  (i=0;  i  <  nno<iesx;  i*») 

( 

for  (j=0;  j  nniMesy; 

f  scan  1  (rnixiif  i  le,  "  i  If  ,  If  ,  ",  imreal,  tmimag); 

mag  ®  fsgrt  (mreal  ♦  mreal  ♦  mimag  ♦  mimag) ; 

phase  “  fatan.’ (mima'i,  mreal); 

phase  '  phase/cdr; 

if  (mag  >  mommax)  m.-mmax  -  ma-j; 

if  (mag  <  mr>mminf  m  mmin  mag; 

mommag [ i | ( i 1  »  maq; 

while  (pliase  <  n.nj  phase  --  phase  t  thO; 

while  (phase  '  ^60.*^)  phase  «  phase  -  360; 

m'-mphase  [  i  )  1  1  )  -  phaf:*^; 

) 

I 

for  (1  =  0;  i  t  nri'xie.ax;  i  +  fi 
I 

for  (]=0;  J  •  r.rio<1**3Y;  i**) 
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d«ta ( i I  I  j )  1 0]  »  2 . 0* i/ (nnodeax- 1 )  -  1.0; 
d«t« [il  I  j 1 ( 1 )  •  2.0* j/ (nnodasy-l)  -  1.0; 

data  ( i  I  ( j  1  (21  <■  (f loglO  (noanin/moaaBaqt  IH3 1 ) /t LoqlO  (mooain/iacMBkaxn  ; 
I 
I 

fclosa(ncinfile) ; 

xmax  •  gatqdeac (GDXPMAXl ; 
ynax  ■  getgdasc (GD_YPMAX) ; 

2val  -  qatgdasc {GD_ZMAX) ; 

minsize (xmax/ 10,  ymax/10); 
max9iza(xiBax-200,  /roax-160); 
kaepaspect (xmax,  /max); 
ao*n_win  •  winopen  (’*••) ; 
wintitle(mom_win_nam«) ; 

winconstrainta ( ) ; 
winaat  (motn  win)  ; 
zbuffer(TRUB) ; 
doublebuffer ( } ; 

RGBmode ( ) ; 
gconf ig ( ) ; 

ahademodeKGOUBAUD) ; 
mmode(KVI  EWING)  ; 

qdevice{LEFTMOUSE) ; 
qdevice(MIDOLEMOUSE) ; 
qdevice(RIGHTMOUSE) ; 
qdevice(ESCKEY) ; 
tielLEFTMOUSE,  MOUSEX,  M'XJSEY); 


1  end  initialized  ♦/ 


/•  main 

*  main  calls  initialize,  then  goea  into  a  loop.  Within  the  loop 

*  drawacene  ia  called  and  eventa  in  the  event  queue  are  processed 

*  until  user  exits  (by  pressing  ESCAPE  or  through  the  window  manager). 

*  Mouse  input  is  passed  to  dcawscene  when  the  left  mouse  button 

*  is  held  down. 


void 

fflain(void) 

( 

Boolean  exitflag  •  FALSE; 

Boolean  redraw^needed  »  TRUE; 

Boolean  draw_moin  ■  TRUE; 

Boolean  dcaw_ph  *  FALSE; 

Boolean  draw_^data  ■  FALSE; 

short  value; 
abort  buffer (BUFSIZEI ; 
long  row__hits  *  0; 
long  col_hits  ■  0; 
short  deltax,  delta/; 
short  row,  column; 

long  dev; 
long  choice; 

long  data  choice,  trans_choice,  display^mode,  iblink; 
long  menu,  datamenu,  tiansmenu,  displaymenu; 

short  X,  oldx;  /♦  mouse  movement  */ 

short  y,  oldy; 

short  dxr,  dxt,  dxz,  dxs  -  0; 
short  dyr,  dyt,  dyz,  dys  »  0; 
static  long  whitecoll)  -  I  255,  255,  255  1; 
static  long  bluecolll  -  (  0,  0,  255  1; 
static  long  dbluecol 1 1  *  I  0,  0,  128  J; 

initialize ( ) ; 
data_choice  *  1; 
trans_choice  •  1; 
display_mode  •  2; 
iblink  >  0; 
tx  »  0; 
ty  •  0; 
tz  •  0; 

/*  make  the  data  menu  */ 

datamenu  »  defpup( '’Choose  Data  Window  to  Open/Close  ttt" 

"Magnitude  txl /Magnitude  «  Phase  1x2"); 

/*  make  the  transformation  menu  */ 
transmenu  •  del  pup ("Choose  Transformation  It  I" 

"Rotate  1x3 /Translate  IxlIZoom  lx5lScale  2-corap  txBlReturn  to  initial  state  IxT") 
/*  make  the  display  mode  menu  */ 
displaymenu  -  defpup("Choose  Display  Mode  Itl" 

"OD  lx9llD  lxlOI2D  txll"); 

/*  make  the  main  menu  */ 

menu  «  defpup("Menu  ItIData  Window  ImtTransformation  ImIDisplay  Mode  Iral" 

"Close  Current  Window  1x8", datamenu, transmenu, displaymenu) ; 


FALSE) 


while  (exitflaq  *** 

1 

if  (redraw  needed  -=  TRUE) 


if  (dfaw_moni) 

drawscene (dxr ,  dxt,  dxz, 

1,  tran3_chc>ice, 
moni_wi  n) ; 

if  (draw_ph} 

drawscene (dxr ,  dxt,  dxz, 

2,  t rans_ch":>ice, 
ph_win/ ; 


dxs,  dyr,  dyt,  dyz,  dys, 
di3play_mode,  ibiink. 


dxs.  dyr.  dyt,  dyz.  dys. 
displaymode,  ibiink. 


if  (ibiink  ==  0)  tedraw_needed  =  FALSE; 
ibiink  =  -ibiink; 

} 


while 

( 


(exit  flag  FALSE)  44 

(qtestO  II  (redraw_needed  =*  FALSE))  ) 

dev  =  qread  (4value); 
if  (dev  ==  ESCKE?) 

{ 

if  (value  ==  0) 


exit  flag  =  TRUE; 

{ 

else  if  (dev  ==  REDRAW) 

( 

reshapeviewport ( ) ; 
redraw  needed  *  TRUE; 

) 

else  if  (dev  »»  INPUTCHANGE) 
i 

our  win  »  value; 

1 

else  if  (dev  RIGHTMOUSE  44  value  1) 

( 

choice  =  dopup<menu); 
if  (choice  >»  1) 

I 

It  (choice  •»  1) 

1 

draw_mom  *  ldraw_mom; 
if  (draw  mom) 

I 

minsize (xmax/ 10,  ymax/lO); 
max3ize(xmax-200,  ymax-160); 
keepaspect (xmax,  ymax); 
mom_win  «  winopen("**)  ; 
wintitle(mom_win_name) ; 

winconstraintsT) ; 
winset(mom  win); 
zbufferfTRUE); 
doublebuf  fer ( ) ; 

RGBmode ( ) ; 
gconf ig () ; 

shaderoodel (GOURAUD) ; 
mniode(MVlEWINC>/ 

1 

else  winclose(mom  win); 

I 

if  (choice  **  2) 


draw_ph  =  !'iraw_ph; 
if  (draw  ph) 

I 

minsire{xmax/10,  ymax/10); 
maxsize (xmax-200,  yraax-lEO); 
keepasp-ect  (xmax,  ymax)  ; 
ph_win  *  winopen("'*) ; 
wintitle(ph_win_name) ; 

winconsttaints () ; 
winset (ph_win) ; 
zbuf fer (TRUE) ; 
doublebuf  fer () ; 

RGBmode ( ) ; 
gconfigO  ; 

shademodel  (CJOURAUD)  ; 
mmode(MVI  EWING)  ; 

1 

else  wxnclo5e(ph_win) ; 


if  {(choice  >=  3)44(choice  <=  7)) 
trans  choice  =•  choice  -  2; 
if  (trans  choice  ==  S) 


dxr  =  0. 
dxt  -  0, 
dxz  •=  0, 
dxs  =  0, 
dyr  «  0, 
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dyt  ■  0; 
dyz  -  0; 
dys  -  0; 
tx  •  0; 
ty  •  0; 
tz  •  0; 

trans_choice  •  1; 

1 

if  (choice  «■  8) 

( 

if  (cur_win  *»  ffion_win)  draw_noei  -  •draw_mcm; 
if  (cur_win  ph_win)  draw_ph  ■  idraw__ph; 
winclo9e(cur_win) ; 


if  {(choice  >-  9)4*(choice  <■  ID)  display_Bode  ■  choice 
redraw  needed  -  TRUE; 


else  if  (dev 


LBFTMOUSE) 


if  (value) 

/*  read  initial  mouse  x,y,  which  will  be 
*  next  2  events  because  tied  to  leftmouse 


{ 


) 

else 

/*  as  soon  as  mouse  released, 

♦  stop  listening  to  mouse  x,y 

( 

unqdevlce  (H(XJSEX) ; 
unqdevice (MOUSEY) ; 


qread (4x) ; 
qread(4y) ; 

/*  now  listen  to  all  mouse  x,y 
*  as  long  as  mouse  down 
qdevice(MOUSEX); 
qdevice (MOUSEY); 


I 

else  if  (dev  ««  MIODLEMOUSE  «4  value  «•  0) 

{ 

/*  do  a  pick  correlation  in  current  window  */ 
winset (cur_win) ; 

/♦  first  find  row  •/ 
pick(buffe£,8UFSI2E); 
picking^row  ■  TRUE: 
deltax  «  3; 
deltay  ■  3; 

picksize(deltax,deltay) ; 
iblink  »  0; 

drawscene(dxr,dxt,dxz,dxs,dyc,dyt,dyz,dy5,data__choice,  trans_choice 
0, iblink, curjwin) ; 
picking_row  -  FALSE; 
row_hits  »  endpick (buffer) ; 
if  (ro/_hits  >  0)  row  «  bufferll); 


/•  now  find  column  ♦/ 
pick(buffer.BUFSlZE) ; 
picking_col  «  TRUE; 
deltax  *  3; 
deltay  «  3; 

picksize (deltax, deltay) ; 
iblink  «  0; 

drawscene(dxr,dxt,dxz,dxs,dyr,dyt,dyz,dy5,data_choice, trans_choice 
0, iblink, cur_win) ;  ' 

picking_col  -  FALSE; 
col_hits  -  endpick(buffer) ; 
if  (col^hits  >  0)  column  -  bufferll); 


if  ((row  hits  >  0) 44 (col_hits  >  0)) 

{ 

if  (draw_data  —  FALSE) 

1 

prefposition(xmax/40,xmax*9/20,yBax/40,ymax*€/40); 

data_win  •  winopen(”") ; 

wintltle (data_win_name) ; 

winconst raints ( ) ; 

winset (data^win) ; 

singlebuffer () ; 

RGBmode  ( ) ; 
gconfigO  ; 
shademodel (FLAT) ; 
maodelHVI  EWING) ; 
draw_data  •  TRUE; 

) 

winset (data_wln) ; 

c3i (dbluecol) ; 

clear () ; 

c3i (whitecol) ; 

ortho2(0.0,  5,0,  0.0,  4.0); 
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printhits (ibiink,  row,  column); 

) 

else 

if  (draw_data) 

1 

draw_data  =  FALSE; 
winclose (data_win) ; 

} 

1 


else 


if 

if 

if 

if 


if  (d^v  **  MOUSEX) 

oldx  “  x; 

X  =  value; 

(tran3_choice  =*  1)  dxr  •  dxr 
{trans_choice  =*  2)  dxt  =  dxt 
(trans_choice  ==  3)  dxz  =  dxz 
{tran3_choice  =*  4)  dxs  =  dxs 
redraw  needed  =  TRUE; 


(x-oldx) ; 
(x-oldx) ; 
(x-oidx) ; 
(x-oldx) ; 


else  if  (dev  ==  MOUSEY) 

I 

oldy  =  y; 
y  =  value; 

if  (tran3_choice  1)  dyr  =  dyr  «  (y-oldy); 

if  (trans_choice  =«  2)  dyt  =  dyt  ♦  (y-oldy); 

if  (t ran3_choice  ==  3)  dyz  *  dyz  ♦  (y-oldy); 

if  (transchoice  4)  dys  »  dya  +  (y-oldy); 

redraw_needed  =  TRUE; 

> 

I  /*  end  while  exittlag  FALSE 

and  (qtest  or  redraw_needed  FALSE)  */ 
)  /•  end  while  exitflag  FALSE  */ 

exit (0) ; 

)  /*  end  mainO  */ 


/*  printhits 

*  Prints  row,  column  info  in  data  window 

♦/ 

void  printhits (long  iblink,  short  row,  short  column) 

I 

int  i , j ; 
char  str{601; 
i  *  row; 
j  »  column; 
cffiov 10.2,3.4,0.0); 

sprintf(str,  ’’This  matrix  has  lid  rows  and  lid  columns.*',  nnodesx,  nnodesy); 
charstr (str) ; 
cmov (0 . 2, 2. 8 , 0. 0) ; 

sprintf(str,  *'You  have  picked  row  41d,  column  lid.",  rcw,  column); 
charstr (str) ; 
cmov (0. 2, 2 .2,0.0); 

sprintf(str,  "The  magnitude  of  the  matrix  at  this  point  is:  4.2f.",  mommag (i ) ( j ) ) ; 
charstr  (str) ; 
cmov (0.2, 1.6, 0.0); 

sprintf(str,  "The  phase  of  the  matrix  at  this  point  is:  t.2f  degrees.",  momphase [ i ) I j ) ) 
charstr  (str)  ; 
cmov (0.2, 1.0, 0.0); 

sprintf(str,  "The  maximum  magnitude  in  this  matrix  is;  mcrttimax); 

charstr (str) ; 
cmov (0.2, 0.4, 0.0) ; 

3printf{str,  "The  minimum  magnitude  in  this  matrix  is:  roommin); 

charstr  (st r  J ; 


/*  h3v2rgb  converts  from  hsv  to  rgb  colors.  It  assumes  that 

*  s  and  V  both  equal  1.  This  can  be  easily  changed  in  the 

*  future.  Result  is  put  into  global  variable  rgbcoll3). 

•/ 

void  hsv2rgb( float  hue) 

I 

float  h, s , V, f , p, q, t ; 
int  j; 

s  *  1.0; 

V  «  1.0; 
h  »  hue/60; 
j  -  h; 
f  =  h-j; 

p  »  V* { 1 -s)  ; 
q  -  V* (1- (s*  f ) )  ; 
t  =  V (l-(s*  (1-f)  ) )  ; 

switch  (j) 

( 

case  0: 

rgbcol (0 ] »v*2S5; 
rgbcol ( 1 ] =t*  255; 
rgbcol 12) =p*255; 
break; 
case  1: 

rgbcol(01=q*255; 
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rgbcol ( I ] *v*255; 
rgbcol (2)“p*255; 
break; 
case  2: 

rgbcol  I0)=p*2S>5; 
rgbcol ( 1 ) ■v*255; 
rgbcol(2]=t*255; 
break; 
case  3: 

rgbcol fO} «p*255; 
rgbcol ( 1 1 »q*255; 
rgbcol [2) »v*255; 
break; 
case  4: 

rgbcol(0)-t*255; 
rgbcol [ 1 ] =p*2SS; 
rgbcol  ( 2 }  **v*  255 ; 
break; 
case  5: 

rgbcol (0) »v*255; 
rgbcol ( 1 ) =p*255; 
rgbcol (2) -q*255; 
break; 


/*  drawscene 

*  drawscene  performs  all  drawing.  Scene  is  cleared  and 

*  entire  scene  is  redrawn  each  time  called-  Mouse  input 

*  is  used  to  change  polarview  viewing  transformation. 


void 

drawscene (short  dxr,  short  dxt,  short  dxz,  short  dxs, 
short  dye,  short  dyt,  short  dyz,  short  dys, 

long  data__choice,  long  trans__choice,  long  display__mode,  long  iblink, 
long  do__win) 

{ 

short  azimuth,  incidence,  fovy; 
float  edr,  azi,  inc,  hue,  xstep; 
float  valmax,  valmin; 
int  i; 

unsigned  precision; 

static  long  v^itecold  •  \  255,  255,  255  I; 
static  long  blackcolU  «  (  0,  0,  0  ); 
static  long  ceded ()  «  |  255,  0,  0  ); 
static  long  orangecoK)  •  |  255,  63,  0  J; 

static  long  yeliowcoljj  -  (  255,  255,  0  J; 

static  long  greencoKI  »  {  0,  255,  0 
static  long  cyancoUl  *  {  0,  255,  255  ); 
static  long  bluecol ( )  »  (  0,  0,  255  1; 
static  long  newbluecoi(J  »  (  0,  52,  255  1; 
static  long  vioietcold  «  (  255,  0,  255  1/ 

static  long  greycolf]  »  j  192,  192,  192  I; 

winset (do_win) ; 
c3i  (blacked) ; 

if  ( (picking^row  ==  FALSE)  <picking_cd  *«  FALSE))  clearO; 

edr  »  M_PI/1800.0; 
fovy  »  450  ♦  dyz; 
if  (fovy  <  10) 

fovy  »  10; 

else  if  (fovy  >  1790) 
fovy  »  1790; 

azimuth  »  dxi  •  5; 
incidence  ■  dyr  ♦  5; 

if  (trans_choice  2) 

( 

azi  >  cdr*azimuth; 
inc  "  edr* incidence; 

tx  «  (dxt* fcos (azi)  -  dyt*fcoa (inc) * fsin{azi) )/100; 
ty  -  (dxt*fsin(azi)  ♦  dyt*fcos(inc)*fco3(azi) )/l00; 
tz  •  (dyt*fsin(inc) )/100; 

! 

perspective (fovy,  5. 0/4.0,  1.0,  9.0); 

if  ( (picking_row  ««  FALSE)  ««  (picking_col  FALSE))  czclear(0x000000,  zval); 
pushmatrix( ) ;  /*  save  ModelView  matrix  by  pushing  and  duplicating  */ 
pdarview(5.0,  azimuth,  incidence,  0); 

/•  viewing  commands  prenultiply  the  ModelView  matrix  */ 
c3i (whitecol) ; 
pushnatrixO ; 

translate(tx,ty,tz) ; 

mommatrix (data_choice,  display^mode,  ibiink) ; 
popmatrixO ;  ” 

c3i (greycol) ; 

if  ( (pickingrow  --  FALSE)  t(  Ipicklnqcol  —  FALSE))  »K«5(); 
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popmat r ix I } ;  /•  restore  ModeiView  mdtrix  by  popping  •/ 

/*  change  projection  to  ortho2  to  display  text  in  2D 
V 

ortho2(0.0,  S.O.  0  0,  4.0); 

it  { {picking_row  ==  FALSE)  &&  (picking_col  ==  FALSE))  swapbuffersi) 
)  /♦  end  drawsceneO  */ 


/***• 

AXES 


/*  three  dimensional  vector  •/ 
typedef  float  vector [3]; 


/*  useful  vectors  for  the  axes  */ 


static 

vector 

front 

- 

{  0.0, 

0.0, 

1.0) 

St  at i: 

vector 

L>ack 

= 

I  0.0, 

u.  0, 

-i.O) 

static 

vector 

top 

1  0.0, 

1.0, 

0.0) 

static 

vector 

bottom 

= 

{  0.0. 

-1.0, 

0.01 

static 

vector 

right 

* 

1  1.0, 

0.0, 

0.0) 

static 

vector 

left 

= 

1-1.0, 

0.0, 

0.0) 

static 

vector 

center 

I  0.0, 

0.0, 

0.0} 

/*  draws  a  wireframe  set  of  right-handed  axes  */ 
void  axes (void) 

( 

bgnl ine ( ) ; 

v3f  (center) ; 
v3f (right) ; 

endline ( ) ; 

cmov (right ( 0] , r  Lght 1 1 1 , r ight ) ) ; 
charstr ("X") ; 

bgnlineO  ; 

v3f (center ) ; 
v3f  (top) ; 

endlineO  ; 

cmov (topic] , topi  1 1 , topi  21 ) ; 
charstr  ("Y") ; 

bgnl ine ( } ; 

v3  f (center ) : 
v3f (front ) ; 

endlineO ; 

cmov (front (01 , front ( I j , front (2] ) ; 
charstr (“2") ; 

) 


MoM  MATRIX 


/ ♦  MoM  matr ix  * / 

void  mommatr ix ( long  data_choice,  long  display_mode,  long  iblink) 
i 

int  i,  j; 

short  rcw,  column; 

float  h, 3 , V, f , p, q, t , vai ; 

static  float  origin!)  =  (  0,  0,  0  ); 

static  long  redcol ( I  =  (  2S5,  0,  0  ); 

static  long  orangecoll)  =  {  63,  0  ); 

static  long  yellowed  []  =  i  266,  2SS,  0  ) ; 

static  long  greened [ ]  =  (  0,  2Sb,  0  1; 

static  long  cysncoll]  -  I  0,  265,  255  ); 

static  long  bluecolf]  =  i  0,  0,  255  ); 

static  long  newbluecol [ ]  =  I  0,  52,  255  ); 

static  long  violetcolM  =  1  255,  0,  255  ); 

static  long  greyed ( J  =  (  192,  192,  192  }; 

static  long  purpled  (1  =  I  255,  0,  176  I; 

static  long  whitecoll)  -  {  255,  255,  255  I; 

if  ( (picking_row  =*  FALSE)  44  \picking_col  =■=  FALSE)) 

I 

if  (display_mode  ==  0) 

I 

bgnpoint ( ) ; 

for  (i»0;  i  <  nnodesx;  i-**) 

( 

for  (3=0;  j  <  nnodesy;  j*+) 

I 

if  (data_choice  l)  hsv2rgb(d3ta  I  i  )(  j  ]  (2]  •  360. 0)  ; 
if  (data_choi‘:e  ^  2)  hsv2rgb(momphase  ( i )  ( j  J )  ; 

c3i (rgbcol) ; 
v3  f (data  1 i ) ( j 1 ) ; 

I 

I 

endpoint ( ) ; 

i 

if  (display_mode  ==  1) 

f 

for  (i»0;  i  <'  nnrdesx;  i**-) 


/ 
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1 

b^nline  ( ) : 

for  (j*0;  j  <  nnodeay;  j++) 

if  (data_choice  ==  1)  h3v2rgb(data [i} [j 1(21 *360. 0) ; 
if  (data_choice  =»=  2)  h3v2rgb(momphase(i)  ( j  J }  ; 
c3i (rgbcol) ; 
v3f (data [i 1  I  j 1 ) ; 

) 

endline () ; 

1 

for  (j=0;  j  <  nnodesy;  j++) 

( 

bgnline ( ) ; 

for  (i*0;  i  <  nnodesx; 

if  (data_choice  =»  1)  hsv2rgb(data ( i It j 1  (21  * 360 , 0) ; 
if  (data_choice  ==  2)  h3v2rgb (momphaaef i J (j I ) ; 
c3i (rgbcol) ; 
v3f  (data ( i] ( j 1 ) ; 

) 

endline ( ) : 

I 

I 

if  (display_mode  **  2) 

{ 

for  (i»0;  i  <  nnodesx-l;  i+*-) 

( 

bgnqstripl ) ; 

for  (5*0;  j  <  nnodeay;  5++> 

I 

if  (data  choice  ==  1)  hsv2rqb(data [ij (j 1(21 *360.0) ; 
if  (data”choice  *=  2)  h3v2rgb{raomphase(il ( jl ) ; 
c3i (rgbcol) ; 
v3f (data (il  [ j ) ) ; 

if  (data_choice  1)  hsv2rgb(data ( i+l 1 ( j 1(21 *360. 0) 
if  (data_choice  =*  2)  hsv2r jbfmoophasel i+ll ( j ) ) ; 
c3i (rgbcol) ; 
v3f (data(i+ll (jl); 

) 

endgstripO ; 

} 

1 

) 

if  (picking_cow) 

I 

for  (i«0;  i  <  nnodesx;  i**) 

1 

loadname(l+l) ; 
bgnpoint {) ; 

for  (j*0;  j  <  nnodesy;  j++) 

( 

v3f(data(inj!); 

1 

endpoint ( ) ; 

) 

) 

if  (picking  col) 

( 

for  (5*0;  j  <  nnodesy;  j*+) 

I 

loadname( j+1) ; 
bgnpoint () ; 

for  (i=0;  i  <  nnodesx;  i++} 

\ 

v3f (data ( i; (51); 

I 

endpoint ( ) ; 

1 

1 
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